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Abstract. In this work we prove global stability for the Boltzmann equation 
(1872) with the physical collision kernels derived by Maxwell in 1866 for the 
full range of inverse power intermolecular potentials, r^^P~^^ with p > 2. This 
completes the work which we began in |33l (and announced in 1341 ). We more 
generally cover collision kernels with parameters s G (0, 1) and 7 satisfying 
7 > — (n — 2) — 2s in arbitrary dimensions T" X R" with n > 2. Moreover, we 
prove rapid convergence as predicted by the Boltzmann H-Theorem. When 
7 + 2s > 0, we have exponential time decay to the Maxwellian equilibrium 
states. When 7+2s < 0, our solutions decay polynomially fast in time with any 
rate. These results are constructive. Additionally, we prove sharp constructive 
upper and lower bounds for the linearized collision operator in terms of a 
geometric fractional Sobolev norm; we thus observe that a spectral gap exists 
only when 7 + 2s > 0, as conjectured in Mouhot-Strain 1491 . 
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1. Introduction, main results, and overall strategy 

The present article considers the Boltzmann equation [T^ from 1872 with the 
collision kernels derived by Maxwell [23 in 1866 for the intermolecular repulsive 
potential 0(r) — r~(^'~^) over the full range p S (2, 00). It will be observed that 
this fundamental equation, derived by both Boltzmann and Maxwell, grants a basic 
example where a range of geometric fractional derivatives occur in a physical model 
of the natural world. 
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This work contributes to the understanding of global in time, close to Maxwellian 
equilibrium solutions of the Boltzmann equation without angular cut-off, that is, for 
long-range interactions. This problem has been the subject of intense investigations 
for some time now. We develop a mathematical framework for these solutions for 
all of the collision kernels derived from the intermolecular potentials and more gen- 
erally. For the hard-sphere and cut-off collision kernels, such a framework has been 
well established for a long time [36H38t.46,,55..56] . The hard-sphere kernel applies in 
the limit when p — >■ oo. A framework is also known for the Landau collision operator 
[55] . which can be thought of as the limiting model for p = 2. The important case 
of the Boltzmann equation with the collision kernels derived by Maxwell is then the 
last case for the intermolecular potentials in which this framework has remained 
open. Our concurrent article |33| established such a framework for the easier "hard 
potential" cases when p > 3. The current article provides this framework for all of 
the intermolecular potentials including some of the most interesting and singular 
"soft potential" cases. We will discuss more of the historical background for the 
Boltzmann equation for long-range interactions in Section 11.21 

The model which is the focus of this research is the Boltzmann equation 

dF 

(1) —+vV^F=Q{F,F), 

where the unknown F{t^ x, v) is a nonnegative function. For each time t > 0, 
F{t,-,-) represents the density of particles in phase space, and is often called the 
empirical measure. The spatial coordinates are a; G T", and the velocities are 
V £ R" with n > 2. The Boltzmann collision operator Q is a bilinear operator 
which acts only on the velocity variables v and is local in (t, x) as 

Q(F, F){v) = [ dv, [ da B{v - w„ a) [F^F' - F^F\ . 

Here we are using the standard shorthand F = F{v), F^, = F(wh.), F' = F{v'), 
Fl = F{v'^). In this expression, v, v^, and u', v'^ are the velocities of a pair of 
particles before and after collision. They are connected through the formulas 

, , , v + v* \v-vA , v + v* \v-v*\ „_l 

(2) V = + ^ cr, = ^ cr, a e § . 

This representation of Q results from making a choice for the parameterization of 
the set of solutions of the physical law of elastic collisions: 

w -f w» — v' + vl, 

\v\^ + K\' = \vT + K\'- 

This choice is not unique, and we also use Carleman-type representations later on. 

The Boltzmann collision kernel B{v~v^, a) for a monatomic gas is a non- negative 
function which only depends on the relative velocity \v — v^\ and on the deviation 
angle 6 through cos9 = {k,a) where k = {v — — u*| and (•, •) is the usual 

scalar product in R". Without loss of generality, we may assume that B{v — v^.,a) is 
supported on {k,a) > 0, i.e., < < ^. Otherwise we can reduce to this situation 
with the following "symmetrization" ; 

B{v - u,,(t) = [B{v - u*,cr) -t- B{v - u*,-(t)] l(fc,ff)>o- 

Above and generally, 1a is the usual indicator function of the set A. 
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The Collision Kernel. Our assumptions are as follows: 

• We suppose that B{v — a) takes a product form in its arguments as 

B{v - u,,cr) ^{\v - v,|) b{cos0). 

In general, both b and <!> are non-negative functions. 

• The angular function 1 b{t) is not locally integrablc; for c;, > it satisfies 



(3) ^<«in"-^^6(cos0)<-^, .G(0,1), V0e(O,- 

Some authors use the notation ly = 2s Cz (0, 2), which is equivalent. 
• The kinetic factor z t-)- satisfies for some C$ > 

(4) ^{\v - v,\) ^ C^\v ~ v,\'> , 7>-(n-2)-2s. 

Our main physical motivation is derived from particles interacting according to 
a spherical intermolecular repulsive potential of the form 

(j){r) = r~^P-^\ pe(2,+oo). 

For these potentials, Maxwell l47| in 1866 showed that the kernel B can be com- 
puted. In dimension n ~ 3, B satisfies the conditions above with 7 = (p — 5)/(p — 1) 
and s = l/{p— 1); see for instance [TBlfrTllSD] . Thus the conditions in ^ and dU 
include all of the potentials p > 2 in the physical dimension n = 3. Note further 
that the Boltzmann collision operator is not well defined for p = 2, see [50] . 

We will study the linearization of ([l} around the Maxwellian equilibrium state 



(5) F{t,x,v) = ^l{v) + ^/fl(v)f{t,x,v), 
where without loss of generality for n > 2 we have 

piv) = (2^)-"/2e-l^l'/2. 

We will also suppose without restriction that the mass, momentum, and energy 
conservation laws for the perturbation f{t,x,v) hold for all t > as 

(6) I dx I dv \ V ] y/JI(uj f{t,X,v) =0. 



This condition should be satisfied initially, and then will continue to be satisfied 
for a suitable classical solution. Our main interest is the global stability of unique 
classical solutions to the Boltzmann equation ([T]) for the long-range collision kernels 
and ^ when the initial data are perturbations of a global Maxwellian ([S]). 
Our solution to this problem rests heavily on our introduction into the Boltzmann 
theory of the following weighted geometric fractional Sobolev space: 

N^'-'^^'ifeS'iR") : l/U.,. <oo}, 
where we specify the norm by 

This space includes the weighted space, for £ € M, with norm given by 

1 2 dof 



\f\h = dv {vY \fiv)\^ 
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The weight is {v) = -^/l + The fractional differentiation effects are measured 
using the following non-isotropic metric d{v, v') on the "lifted" paraboloid: 

d{v,v') = ^J\v~v'\^ + ^{\v\^-\v'\^f. 

The inclusion of the quadratic difference | w p — | p is essential; it is not a lower order 
term. Heuristically, this metric encodes the non-isotropic changes in the power of 
the weight, which are entangled with the non-local fractional differentiation effects. 

The space N'^'^ is essentially a weighted non-isotropic Sobolev norm; this partic- 
ular feature was conjectured in [33]. We see precisely that if M" is identified with a 
paraboloid in R"+^ by means of the mapping v H> (w, and Ap is defined to 

be the Laplacian on the paraboloid induced by the Euclidean metric on R"+^ then 

I/Ik. « / dv{vy+''\{I~Ap)if{v)\\ 

This holds in the sense that the sharp geometric Littlewood-Paley characterization 
of this space is the same as that of N''''^ . The rest of our Sobolev spaces are defined 
in Section ll.il iust below. We may now to state our main result as follows: 

Theorem 1. (Main Theorem) Fix K > K[n), the total number of derivatives, 
and ^ > the order of the velocity weight in our Sobolev spaces. Choose initially 
fo{x,v) G H^iT"" X R") in ([5]) which satisfies There is an rjo > such that if 
||/o||_f/« ^ Vo, then there exists a unique global classical solution to the Boltzmann 
equation ([T]), in the form ([5]), which satisfies 

f{t,x,v) e LfHf{[Q,oo) X T" X R") n L2Ar;;^((0, oo) x T" x M"). 

When J + 2s > 0, for a fixed X > 0, we have exponential decay as 

ll/(*)llHf (T"xK") ^ ^"^'ll/oll/i-if (T"xK")- 

When 7 + 2.S < 0, if fo e H^^„^(T" x M") for weights £,m>Q then we have 

ll/(*)llHf (T"xR") <Cm(l + f) '"||./o||h|^^JT"xK")- 

We also have positivity, i.e. F = ii + ^if > if Fq = fi + y^fa > 0. 

For the derivatives used above, in dimension n = 3, we use K{3) ~ 5. Generally 
K{n) is the number of derivatives required by our use of the Sobolev embedding 
theorems in n dimensions. Here K{n) = 3K* — 1 is sufficient, where K* = [f + Ij 
is the smallest integer which is strictly greater than ^ . This regularity is not needed 
for the hard potentials as the theorem below shows: 

Theorem 2. (Main Theorem from |33| j Fix Kx > 3, the total number of spatial 
derivatives, and Ky > the total number of velocity derivatives. Suppose j+2s > 0. 
Choose initially fo{x,v) € H^'-' (M.^; H^^ (T^)) in ^ which satisfies dH). There is 
an rjo > such that if \\,fa\\H'^v(K.^:H'^^(T^}} ^ VOj then there exists a unique global 
strong solution to the Boltzmann equation in the form which satisfies 

f{t,x,v) e Lt°°([0,oo);i?^"i/^-(T3 x R^)) n L^ti{0,oo); N^'J^^H^^^- {T^ x R^)). 

Moreover, we have exponential decay to equilibrium. For some fixed A > 0, 

\\f{t)\\H'<v(M.3.HK^{T3)) < e"^*||/o|l/i-A'„(R3.^K,(T3)). 

We also have positivity, i.e. F = fj, + ^Jfxf > if Fq ^ fi + ^/Jlfo > 0. 
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We do not attempt to optimize the amount of regularity needed in the above 
theorems, ahhough a trick used just after Lemma |9] allows us to reduce the standard 
regularity conditions (as required for instance by Guo [321I3S] for the cut-off soft 
potential and the Landau equation) from /v (3) = 8 down to K{3) = 5, etc. 

Before we discuss these results wc will now define the relevant Sobolcv spaces. 

1.1. Function spaces. Wc first define a unified weight function 

, , dcf / {v) , 7 + 2s > 0, "hard potentials" 
^ " I {vy^'^' , 7 + 2s < 0, "soft potentials". 

This somewhat non-standard terminology distinguishes exactly when a spectral gap 
exists for the linearized collision operator. In both cases the weight goes to infinity 
as |f I goes to infinity. The scaling chosen above will be convenient for the soft 
potentials in our analysis below. 

We will use the multi-indices a ~ (a^, a^, . . . , a") and /3 = (/3\ . . . , 
Inequalities for multi-indices are defined as follows: if each component of a is not 
greater than that of ai, we write a < ai. Also a < ai means a < ai and 
|a| < |q;i|, where |q;| = -f + • • • + a" as usual. We will use the spatial and 
velocity derivatives 

= dtdt ■ ■ ■ dfA\€. ■ ■ ■ dC ■ 
For any ^ G M, the space i7^(R") with K > velocity derivatives is defined by 

\P\<K 

The norm notation |/i|?rK '^\b\<k |w^^'^'i9/3/i|?2 ^jjti-i will also find utility. 

It will be convenient to define the total H ~ iJ/^(R" x T") space by 

5^^lli2(T"xK")- 



\cx\ + \l3\<K 

We also consider the general weighted non-isotropic derivative space as in ([7]) by 

|.|2 dcf I £j.|2 , /" J / \7+2s+l 21/ \ f -I I {f W) ~ f i""))^ -, 

\h\Mr - 1^ /L^,.. + dv {v) w (v) J^^ dv 

The L'^(]R")-based norms exclusively in the velocity variables are always denoted 
as I • I with the space indicated by the subscript. The usual L^(R") inner product is 
then (•,•), which because of the context should not be confused with our notation 
for the scalar product. We will define L^(T" x M") norms in both space and 
velocity variables by replacing | • | with j| • j| and keeping the same norm subscript. 
We also use (•, •) to be the standard L^(T" x K") inner product. In particular, 

II^IIat-.t = II \h\N'n Wl^Tn.) and ||/i|U;-^ = lll'iU;'T-||L2(T")- We wih also use the 
high derivative, non-isotropic space N = Ng'^{T'^ x M"), given with K > hy 

II^IIa = ll^llAr;'^(T"xR") = E ll^/J^llw;i^^|(T"xR")- 

\a\ + \f!\<K 

Furthermore I^Iat^'^ = |^lAf°'^(R") is the same norm as above, but only in the veloc- 
ity variables. We remark that the decaying weight lu"!'^' on the velocity derivatives 
is only needed for the soft potentials. However we may use it for free for the hard 
potentials and it does not reduce the size of the space. Additionally, in Theorem [2] 
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above the space -/V^'^ denotes specifically the regular N"''* space with Ky isotropic 
velocity derivatives and no weights. Also in Theorem [2] the space denotes a the 
usual Sobolev space with K derivatives in either the x variables or the v variables. 

1.2. Historical Discussion. For early developments in the Boltzmann equation 
for long-range interactions, between 1952-1988, we mention the work of Arkeryd, 
Bobylev, Pao, Ukai, and Wang Chang-Uhlenbeck-de Boer in [UTTllSnilSIlini] • 

Otherwise, many mathematical works have required the Grad |32| angular cutoff 
assumption from 1963-6(cos0) G L°°{S'^) or sometimes b{cos9) G L-'^(S^) -which 
is designed to avoid the non-integrable angular singularity from ([3]). We refer the 
reader to a brief few breakthrough works of pl l [T7 l [26 l [27 l [3CT l l32 l[MgTl l43 l l46 l l52l^ ; 
further references can be found in the review article |60| and |33| . Notice that many 
of these works develop ideas and methods that are fundamental and important 
even without angular cut-off. In particular the space-time estimates and general 
non-linear energy method developed by Guo [351 - 137] and Strain-Quo [53l[54] is an 
important element in the last section of our proof. 

Yet with the exception of the hard sphere model, which is derived in the limit 
when p — > oo, the rest of the inverse power-law potentials dictate that the cross 
section B{v — v^,, a) is a non-integrable function in the angular variable and the 
collision operator, Q{f,f), is then a non- linear diffusive operator. 

In particular, in the presence of the physical effects due to the angular sin- 
gularity, the Boltzmann equation is well-known to experience regularizing effects. 
These results go back to Lions [31] and Desvillettes [TO] and have seen substan- 
tial developments [2l [20H22l[24l[25l[59] . Recently Chen-Desvillettes-He [18] and also 
Alexandre-Morimoto-Ukai-Xu-Yang [6l[7] have developed independent machinery 
to study these general smoothing effects for kinetic equations. 

Contrast this with the case of an angular cut-off, where, as a result of works by 
Boudin-Desvillettes [2] in 2000, and additional progress in |T^[55], we know that 
under the angular cut-off assumption small-data solutions have the same Sobolev 
space regularity as the initial data. The results described in these last two para- 
graphs illustrate that the Boltzmann equation with angular cut-off can be in some 
respects a very different model from the one without any angular cut-off. 

We also mention several works which developed and utilized the entropy pro- 
duction (or dissipation) estimates for the collision operator (which grants a non- 
linear smoothing effect) and also the spatially homogeneous theories. We point out 
the work of Alexandre. Bouchut, Desvillettes, El Safadi, Fournier, Golse, Goudon, 
Guerin, Lions, Mouhot, Villani and Wennberg as in [2l|4l[5] [22H24l l29 p 3i p 45 p 58j . Fur- 
ther references can be found in our concurrent article [33] and the surveys [3ll60]. 

Lastly, for the most physically interesting and mathematically challenging case 
of the spatially inhomogeneous Boltzmann equation there are much fewer results. 
Here we have two results on local existence [HIZ]. We also refer to the paper by 
Alexandre- Villani [5] , which proves global existence of DiPerna-Lions renormalized 
weak solutions [27] if one includes in the equation a non-negative defect measure. 

Furthermore, after the completion of this present work and |33[I34| . we were 
kindly informed by R. Alexandre of the recent paper [5] by Alexandre-Morimoto- 
Ukai-Xu-Yang. They prove global existence and smoothness for perturbations of the 
Maxwellian equilibrium states ([5]) in the whole space R'^ for the Maxwell molecules 
collision kernel (meaning that the kinetic factor in (|4]) is constant) and moderate 
angular singularities (meaning that 0<s<l/2in([3])); these assumptions apply 
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to the inverse power intermolecular potentials when p = 5. The methods of proof 
in [8] are substantially different from our own. 

For the Boltzmann collision operator, its essential behavior has been widely 
conjectured be the same as a fractional flat diffusion — (— A„)*. Precisely 

F ^ Q{g, F) - -(-A^)-'F + l.o.t. 

Above "l.o.t." indicates that the remaining terms will be lower order. The original 
mathematical intuition for this conjecture has been credited to Carlo Cercignani 
[TS] in f969, now more than forty years ago (see for instance Villani [60l p. 91]). 
This conjecture is certainly correct in terms of the smoothing effects, as has been 
proven for instance in [linillKIlHlllinilllini . 

Even so, our research proves at the linearized level that the essential behavior of 
the collision operator is not exactly a flat fractional diffusion. Instead, the fractional 
diffusive effects are fundamentally intertwined with non-isotropic velocity weights 
which go to infinity at infinity. We prove that this interconnection can be charac- 
terized geometrically; the sharp linear behavior is in fact that of a fundamentally 
non-isotropic fractional geometric Laplacian, the geometry being given by that of 
a "lifted" paraboloid in The intuition for this behavior is derived from the 

original physics representations for the collision operator in terms of (5-functions. 

Our next estimate below further establishes that this fundamental non-isotropy, 
with the same geometry, carries over to the fully non-linear situation in many cases. 

1.3. Entropy Production. Among Boltzmann's most important contributions to 
statistical physics was his celebrated i7-Theorem. We define the H-functional by 



H{t) = - [ dx j dv F log F. 

idicts that the entri 
[ dx D{F) > 0, 



Then the Boltzmann H- Theorem predicts that the entropy is increasing over time 

dH{t) 
dt 

which is a manifestation of the second law of thermodynamics. Here the entropy 
production functional, which is non-negative, is defined by 

D[F)^-( dvQ{F,F)\ogF 

JR" 

= - [ dv [ dvJ d<jB{v - V,, a) (F'Fi-FF,) log 

Moreover, the entropy production functional is zero if and only if it is operating on 
a Maxwellian equilibrium. These formulas formally demonstrate that Boltzmann's 
equation defines an irreversible dynamics and predicts convergence to Maxwellian 
in large time. Of course these predictions are usually non-rigorous because the 
regularity required to perform the above formal calculations is unknown at the 
moment to be propagated by solutions of the Boltzmann equation in general. 

In particular, even though there are many important breakthroughs in this di- 
rection, none, so far as we are aware, can be said to completely and rigorously 
justify the H-Theorem for the inverse power-law intermolecular potentials with 
p € (2,oo). We refer to ISl[ll[ll[171[Sl[Sl[SD] , and the references therein, in this 
regard, to mention only a brief few works. Our results in this paper prove rapid 
convergence to equilibrium for all of the inverse power collision kernels derived by 
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Maxwell when one considers initial data which is close to some global Maxwellian. 
This convergence is the essential prediction of Boltzmann's i7-Theorem. 

Many works study entropy production estimates in the non cut-ofF regime, as in 
[4ll45| . These estimates have found widespread utility. Our results and non- isotropic 
Sobolev space have the following implications in the fully non-linear context. Com- 
bining estimates in |33 1 l59 j we can see the following lower bound 



DiF)>[ dv [ dv' i{v){v'))^ i,(.,.,)<i-i-"-t- 

Precisely, this estimate follows straightforwardly from our new estimates [331 Sec- 
tion 6] when combined for instance with the decomposition |591 (29)] in the partic- 
ular case when the unknown functions in the term A{v,v') from [59l (29)] satisfy 
F(v) > £ e"'^'"' for any given £,c > 0. Remarkably, this is the same semi-norm 
as in the linearized context, and it is a stronger non-isotropic and non-local ver- 
sion of the local smoothing estimate from [4] (stronger, that is, in terms of the 
weight power multiplied on the order of differentiation). We have not as of now 
attempted to prove this estimate away from the regime where F is bounded below 
by a Maxwellian. This estimate was derived as a result of our effort to find explicit 
equivalence between the non-isotropic norm coming out of the linearized collision 
operator and the non-isotropic Sobolev space N'^'^ . 

Our methods introduce a new tool into the study of the Boltzmann equation 
with non-intcgrable angular cross-sections ([3]) which allows to us to obtain sharp 
control over geometric fractional differentiation effects with respect to both the 
order of differentiation and the non-isotropic geometrically dependent weight. We 
expect this approach to be useful for future studies in the Boltzmann theory. 

Acknowledgments. Recently Carlo Ccrcignani has passed away in January of 
2010, we wish to take this opportunity to express our deep gratitude and respect 
for all that he has done to advance the field of kinetic theory. 

1.4. Reformulation. Throughout this sub-section we will define the relevant no- 
tation. Wc also reformulate the problem in terms of the equation ([S]) for the pertur- 
bation ([5]). In the next sub-section, we will state and explain our main estimates. 

Throughout this paper, the notation A < B will mean that a positive constant 
C exists such that A < CB holds uniformly over the range of parameters which 
are present in the inequality (and that the precise magnitude of the constant is 
irrelevant). In particular, whenever either A ov B involves a function space norm, 
it will be implicit that the constant is uniform over all the relevant elements of 
the space unless explicitly stated otherwise. The notation _B > ^ is equivalent to 
A < B, and A ^ B means that both A < B and B < A. Furthermore, for any 
i? > 0; we will use Bj^ C M" to be the ball centered at the origin of radius R. 

Wc linearize the Boltzmann equation ([1} around the perturbation ([5]). This 
grants an equation for the perturbation f{t,x,v) as 

(8) dtf + v-V,f + Lif)^Tif,f), fiO,x,v)^fo{x,v), 

where the linearized Boltzmann operator L is given by 

Lig) - - f^-'^'Qifi, VJ^g) - m"'/'Q(Vm.9, /^) 

dv^ I da B{\v - w*|,cos6l) [g^M + ghh - g'M ~ 9 K] M*, 
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and the bilinear operator F is given by 

(9) r(.g,/i) =V"'/'Q(VM5,V7^/^) = J dv, daBAhig'.h' - g,h). 
In both definitions, we take 

M{v) = v/^R = (27r)-"/4e-l"lV4, 

Finahy, we note that 

(10) L{g)'^'-r{M,g)-r{g,M). 

This formulation shows that it is crucial to obtain sharp estimates for F. 
We expand the main term of the linearized Boltzmann operator as 

T{M,g) = j dv, da B{\v~v,\, cose) [Mig' - M,g]AU. 

We will now split this in parts whilst preserving the cancellations as follows: 

F(M, g)^ f dv, f daB [g' - g) M^M, ~ y{y) 9{v). 

where 

i>{v) = / dv^ da B{\v - cos6') (A4 - Mi)M*- 

The first piece above contains a crucial Hilbert space structure. This can be seen 
from the pre-post collisional change of variables as 



dv I dv^ I da B{\v ~v^l cos 9) {g' - g)hMlM^ 

-\ I dv I dv^ [ da B [g - g)hMlM^ 
-( dv I dv^ I da B{g - g')h'M^Ml 

2 JR" JR" JS"-! 



dv / dv^ / da B {g' - g){h' - h)MlM^. 



2 

^ 2 

For the weight, we will use Pao's splitting as 

v{v) = v{v) + vk{v), 

where under only ([3]) and ([4]) the following asymptotics are known: 

v{v) K {v^^^' and \yK{v)\<{vy. 

These estimates were established by Pao in [50l p. 568 eq. (65), (66)] by reducing 
to the known asymptotic behavior of confiuent hypergeometric functions. 

We further decompose L = N + K . Here N is the "norm part" and K will be 
seen as the "compact part." The norm part is then written as 

(11) Ng = -T{M,g) - yK{v)g = - / dv^ f da B{g' - g)M'Jh + i^{v)g{v). 



From previous calculations, this norm piece satisfies the following identity: 
{Ng,g) = lf dv f dvJ daB{g' - gfM^M, + f dv i^iv) \giv)\' 

^ JR" JR" JR" 
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As a result, in the following we will use the non-isotropic fractional semi- norm 

(12) |<7||/= ^ / dvw\v)( dv, [ d<j B [g' - gfMlA'h, £ S M. 

For the second part of {Ng,g) we recall the norm 1/1^2^^ defined below equa- 
tion (O. These two quantities define our norm, which characterizes precisely the 
linearized Boltzmann H- Theorem. 

We also record here the definition of the "compact piece" K: 

(13) Kg = VK{v)g~T{g,M)^VK{v)g-( dv, [ daBAhig'^M' ~ g^M). 

JR" JS"-1 

This is our main splitting of the linearized operator. We will now state our main 
estimates for each of these decomposed operators. 

1.5. Main estimates. In this sub-section we will state most of the crucial long- 
range estimates to be used in our main results. In the next sub-section, we discuss 
how these estimates in particular resolve a conjecture from Mouhot-Strain |49| . 

We will prove all of our estimates for functions in the Schwartz space, 5(]R"), 
which is the well-known space of real valued C°°(R") functions all of whose deriva- 
tives decay at infinity faster than the reciprocal of any polynomial. Note that the 
Schwartz functions are dense in the non-isotropic space N^'J^, and the proof of this 
fact is easily reduced to the analogous one for Euclidean Sobolev spaces by means 
of the partition of unity as constructed, for example, in our first paper |33| . More- 
over, in all of our estimates, none of the constants that come up will depend on 
the regularity of the functions that we are estimating. Thus using routine den- 
sity arguments, our estimates will apply to any function in N'^^'J. or whatever the 
appropriate function space happens to be for a particular estimate. 

The first important inequality that we explain is for the linear operator: 

Lemma 3. (Linear estimates) Consider the linearized operator L — N + K where 
N is defined in ()11|) and K is defined in ()13|) . We have the uniform inequalities 

(14) \{w^'Ng,g)\<\g\l..-,, 

(15) \{w^'Kg,g)\ < vl^'gll^^^^ + C,\g\h^Sa^), 

where £ € R, rj > is any small number, and C,, > 0. 

In these estimates there are several things to observe. First of all, there are 
no derivatives in the "compact estimate" from ()15p . which can be contrasted with 
the corresponding estimate in the Landau case |35[ Lemma 5], in which the upper 
bound requires the inclusion of derivatives. Further p4)) tells us that the coercive 
inequality for the "norm piece" of the linear term in Lemma |4] is essentially sharp. 

Lemma 4. (Main coercive inequality) For the sharp space defined in ([7]) with (jlip. 
and g]) we have the uniform coercive lower bound estimate: 

{w^'Nf, f) > l/l^,.. - C|/|i.(5^), 3C > 0. 

This holds for any £ Cz M.; if £ = we may take C = 0. 

We turn to our first non-linear estimate: 
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Lemma 5. (Non-linear estimate) Consider the non-linear term ([9]). For any \a\ + 
|/3| < K , with K > Kin) and ^ > 0, we have the following estimate 

w^'-^\^\d$T{g^h),d'^pf)\<{\\gU\\h\\H + \\g\\H\\h\\^ 

Notice that this estimate does not require a velocity weight which goes to infin- 
ity at infinity. This feature is made possible by our non-isotropic Littlewood-Paley 
adapted to the paraboloid, which characterizes exactly the geometric fractional dif- 
ferentiation effects that arc induced by the linearized Boltzmann collision operator. 

Lastly, we have two coercive interpolation inequalities for the linearized operator: 

Lemma 6. (Coercive interpolation inequalities) For any multi-indices a,/3, any 
£ > 0, and any small 77 > there is a positive constant Crj, such that we have the 
following coercive lower bound for the linearized collision operator 

(16) {w''-^\^\d"pLg,d^g)>\d^g\%.^.^^-r^ " a,|a"5li.(Bo,)- 

l/3i|<l/3| 

Furthermore, when no derivatives are present, for some C > we have 

(17) {w^'Lg,g)>\g\l.^.-C\g\l.^Scr 

This concludes our statements of the main estimates that will be used in Section 
[5] to establish our main Theorem [TJ Next, we deduce some consequences of these 
estimates for the spectral properties of the linearized Boltzmann collision operator. 

1.6. Conjecture from Mouhot-Strain [JHI- We will now discuss sharp construc- 
tive coercivity estimates of the linearized collision operator, L, away from its null 
space. More generally, from the H-Theorem L is non-negative and for every fixed 
(t, x) the null space of L is given by the (?? 4- 2)-dimensional space 

(18) M = span {y^, vi^, Vn^/Ji, \v\^ ^/Jj] ■ 

We define the orthogonal projection from L^(R") onto the null space N by P. 
Further expand P/i as a linear combination of the basis in (|18p: 

(19) Vh= \ a'\t,x) + Y,b''j{t,x)v, + c''{t,x)\v\A ^. 



;|{I-P}5|^.. < {Lg,g)<C\{\-J>}g\l. 



It is well-known for the linearized collision operator ()T0|) that Lh = if and only 
if ft, = Pft; see the discussion in [33] and the references therein. By combining the 
constructive upper bound estimates in Lemma [3] with the constructive coercivity 
estimate from Theorem [Tl] in Section [Sj we obtain 
1 

C' 

with a constant C > that can be tracked from the proof. Thus a spectral gap 
exists if and only if 7 + 2s > 0, as it was conjectured in [33] ■ The main tools in our 
proof of these statements are the new constructive estimates from Lemma [3] and 
Lemma m combined with the constructive but non-sharp coercive lower bound from 
[48] for the non-derivative part of the norm. This may be of independent interest. 

Notice that for the linearized Landau collision operator, this statement already 
had been shown several years earlier in [35l|48l|49] for any 7 > —3. For Landau, 
there is a spectral gap if and only if 7 + 2 > 0; the Landau operator can be thought 
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of as the limit case when s = 1 and regular (rather than fractional) derivatives are 
present. 

1.7. The proof of Theorem [2l We will now explain how to extend the results 
of |33j to include the more general assumptions stated in Theorem [2] The main 
theorem of [33] gives our Thcorcm[2]for 7 > — min{2s, 3/2}, s e (0, 1), and Ky = 
in three dimensions. Notice that all of the estimates in [33] hold for any 7 + 2s > 
save for the second cancellation estimate in [331 Proposition 6]. However, following 
the remark in the proof of this proposition, we can easily use ([58)) with i = K = 
to extend this lone proposition to all 7 + 2s > 0. The only additional obstruction is 
that the proof of the coercive estimate in [331 Lemma 6] does not directly generalize 
to 7+2s = 0. However notice that TheoremfTTjherein is both constructive and holds 
in full generality of 7 + 2.s > — (n — 2), we can thus use this new estimate in place 
of [33[ Lemma 6]. Lastly, notice that one can further include velocity derivatives 
Ky > in the proof of the main theorem of [S^ and our Theorem [2] using (|22|) . All 
of the proofs in [33] remain exactly the same when velocity derivatives are included, 
except for the cancellation estimates. However the cancellation estimates from [33] 
can be directly extended to cases with derivatives using the decompositions at the 
beginning of the proofs of Proposition [5] and Proposition [5] herein. The rest of the 
proof of the the cancellation estimates in [33| with velocity derivatives are exactly 
the same. 

1.8. Outline of the article, and overview of our proof. The methods from 
our first paper |33j are the starting point for our proof. In particular, several 
important concepts, such as the space N'^''^, the dual formulation for {T{g,h),f), 
the coercive inequality for the semi-norm ([12]), and the non- isotropic Littlewood- 
Paley theory adapted to the paraboloid, will be crucial to our analysis. However, the 
full range of kinetic singularities from ^ dealt with in this paper present several 
additional difficulties not present in [33] and the resolution of these difficulties 
requires numerous widespread essential changes. For example, the L^(R^) approach 
from [33 fails because the kinetic factors from ([4]) are too singular. Instead, in this 
work we include velocity derivatives to control the strongly singular kinetic factors 
via the Sobolev embedding theorems. The addition of velocity derivatives acting 
on the Boltzmann operator requires several new developments, particularly in our 
cancellation estimates, where we provide a new five-term decomposition of the the 
relevant differences. 

Another new theme, in contrast to [33], involves the inclusion of velocity weights. 
Positive velocity weights are used to obtain the rapid decay as in ^3]. Negative 
weights are used to close our estimates, when 7 -I- 2s < 0, as in [55]. Both of 
these uses were known in the cut-off regime; however, the inclusion of these weights 
require several new developments when angular singularities are present. A crucial 
theme in this respect will be to extract a sufficiently large amount of velocity decay 
from the Boltzmann operator, which should be embedded in the functional spaces 
that form the upper bounds in all of our estimates. This is, of course, related to 
the fact that the norm N''''' from ([7]) contains velocity weights which go to zero at 
infinity when 7 + 2s < 0. 

The rest of the paper is organized as follows. In Section [21 we decompose the 
angular singularity and prove the main estimates for each of the decomposed pieces; 
these include size and support estimates, cancellation estimates, and also a useful 
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"compact estimate" in the spirit of Grad which is new to tlie non cut-off theory. 
This compact estimate is the key new element which yields our resolution of the 
conjecture from |49] . Then in Section |3] we develop the non-isotropic Littlcwood- 
Paley theory adapted to the paraboloid in arbitrary dimensions n > 2. As before, 
we have drawn inspiration from the earlier works of Stein |51| and Klaincrman- 
Rodnianski [42], though we arc able to use a more concrete formulation of our 
Littlcwood-Paley projections than is possible in either of these cases because of the 
explicit nature of the underlying geometry. A key point here, new in comparison 
to our first paper [53], is to develop Littlewood-Paley projections that have higher 
regularity than was needed in [33], which requires a more explicit analysis of the 
necessary moment conditions imposed on the kernels. This is crucial for our high- 
order isotropic velocity derivative estimates of the non-isotropic Littlewood-Paley 
projections. After that, in Section |4] we use the estimates from the previous two 
sections to sum all of the decomposed pieces of the Boltzmann collision operator 
from the previous subsections. We sum over the both decomposed singularity and 
the unknown functions expanded via the geometric Littlewood-Paley. In particular, 
we finish the proof of all of the estimates which are stated in Section 11.51 Finally, 
in Section [5l we use methods from the non-linear energy method as in Guo [37j and 
Strain- Quo [S3] to establish global existence and rapid decay. 



In this section, we introduce the first major decomposition of the collision op- 
erator and prove several estimates which will play a central role in establishing 
many of the main inequalities stated in Section 11.51 This is a decomposition of the 
singularity of the collision kernel. For several reasons, it turns out to be useful to 
decompose ^(cos^) from ([3]) to regions where 6 k, 2~^v — rather than a sim- 

pler dyadic decomposition not involving |u — The primary benefit for doing so 
is that this extra factor makes it easier to prove estimates on the space L^_|_2s(M") 
because the weight — ?;*|) from ^ is already present in the kernel and the 
extra weight \v — v^\'^^ falls out automatically from our decomposition. 

The estimates to be proved fall into two main categories: the first arc various L^- 
and weighted L^-incqualitics which will follow from the size and support conditions 
on our decomposed pieces. The second type of estimate will assume some sort of 
smoothness and obtain better estimates than the earlier estimates by exploiting the 
cancellation structure of the non-linear term F from ([9]) . It is already worth stating 
at this point that the particular smoothness assumptions we make are dictated by 
the problem and will specifically be somewhat unusual; in particular, they will not 
correspond to the usual Euclidean Sobolev spaces on R". 

In the rest of the paper, we will carefully explain various arguments the first 
time they are used and subsequently use them again and again with less detail, but 
with references to their earlier appearances. This will save having rewrite similar 
arguments numerous times. 

2.1. Dyadic decomposition of the singularity. Let {Xfej^L-oo ^ partition 
of unity on (0, oo) such that |x/c|l=^ < 1 and supp(x/£) C [2~'^~^, 2"*^]. For each fc. 



2. Physical decomposition and cancellation estimates 
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Notice that we have the expansion 

With this partition, for any ^ e M and any multi-index j3 we define 

Tl'^g, h, f) = I dv ( dv, I da Bk{v -v,,a) g,h Af^«) /' w^'{v'), 

Tj {g,h,f)^ / dv dvr / da Bk{v ~ v^,,a) g^^h Mp{v^,) f w'^^{v). 

We use the notation M^(w) = S/jM = pj3{v)M{v) where ppiv) is the appropriate 
polynomial of degree The dependence of these operators T_^' on /3 is left 
implicit for notational simplicity; in particular, all of the estimates below hold for 
any multi-index /3. 

It will also be necessary to express the collision operator ^ using a "dual for- 
mulation" which is derived via a Carleman-type representation in our first paper 
[551 Appendix] . We record here the following alternative representation for T^'^ as 
well as the definition of a third trilinear operator T^'^ as follows: 

/oi\ rr.k,i /" J / 2f/ /■ , f , Bk g^f'Mfi{v'^)h 

Jr" Jk" Jb^^; 2^ " - - ^ 

Above T^'^ = T^'^{g, h, /), T^'^ = T^'^{g^ h, /), and we use the notation 

B, = i^.(.-..2.'_.-..) = ((^, g::::::, )) x.(k-'i). 

In these integrals diTv is Lebesgue measure on the {n — l)-dimensional plane 
passing through v' with normal — , and v is the variable of integration. More- 
over, above v'^ = v + — v' and M'^ = on . Notice that the insertion of 
the partition of unity guarantees that the kernel is locally integrable for any k. 

Recall the bilinear collision operator © and the collisional variables ([5]). With 
the change of variables u = — w, and = (it ± |it|cr)/2, we may write © as 

r{g,h)= du daB{\u\,a)M{u + v){g{v + u+)h{v + u^) - g{v + u)h{v)} . 

Differentiating this formula and applying the inverse coordinate change allows us 
to express derivatives of the bilinear collision operator F as 

(22) d^Tig,h)= E C^if^rp,{d^Z^^g,d^lh). 

Pl+p2=fi ai<a 

Here C^'^\'^^ is a non-negative constant which is derived from the Leibniz rule. 
Also, is the bilinear operator with derivatives on the Maxwellian M given by 

r^(g, /i)= / dv^, / da B(\v ~ v^,\,a)Mi3{v^,){glh' ~ g^h). 
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Then for ,f,g,h e 5(M"), the pre-post collisional change of variables, the dual 
representation |331 Appendix], and the previous calculations guarantee that 

oo 

k— — oo 
oo 

= E {Tt-\9,hJ)^Ty{g,hJ)}. 

k— — oo 

These will be the general quantities that we estimate in the following sections. The 
first step is to estimate each of T^'^, T^'^, and using only the constraints on the 
size and support of Bk in Section [221 Then in Section [^751 we estimate the collision 
operator again, this time also exploiting the cancellation properties of F^. Finally, 
in Section 12.41 we prove the "compact estimates" which are important in proving 
the constructive lower bound for the linearized collision operator in Theorem 1111 

2.2. Size and support estimates of the decomposed pieces. We will now 
prove several size and support estimates for the decomposed pieces of the Boltzmann 
collision operator. It will be useful to let 0(f) denote an arbitrary smooth function 
which satisfies for some positive constants C'^ and c that 

(23) |0(f)| < C^e-^l"!'. 

We use generic functions satisfying ([23]) often in what follows. 

Proposition 1. For any integer k, any m > and ^ G M, we have 



(24) 
(25) 



Tt^g^hJ) <2'''\g\L^_Jh\^., 



7li 



Furthermore, for cj) defined as in (j23p . we have 



(26) 



T-''(5, </-,/) + r^^'(5,/,' 



<Cs 2 



2sfe 



\9\l^ 



These estimates hold uniformly. 

We record here the following the Sobolev embedding theorem 



(27) 



< 



I/I,.'. 



which holds for any m G M. Wc will also use the following immediate implications 



\w'f\i 



< 



\w-"^.f\L^ + \w'f\L.<\f\^.^, 



where m is sufRciently large, depending on ^ G K, 7 + 2s, and n > 2. Furthermore, 
we are using the abbreviated notation L°° = L°°(IR"). 

Proof. Given the size estimates for b{cos9) in Q and the support of Xk, clearly 



(28) 



da Bk < $(|w - 



de 0- 



-l-2s 



< 2 



\v — v^ I 



This also uses (g]). Thus 



T^j\g.h,f) 



dv 



dv, v/^K'-«*r+^l5*^/l w'^^iv). 
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To obtain ((24|). we first take the L°° norm of g^M^'^ and use (|27l). then apply the 

elementary inequality dv„ mI^'^ |u — u*]''''*'^''* < (v)'^'^'^'^ , and lastly use Cauchy- 
Schwartz putting h in one term, / in the other, and half of the rest in both. 

For ([26| notice that if either of the second two functions in T^'^ has rapid decay 



such as ([23|l . then we have rapid decay in both v and simultaneously. Thus we 
can use Cauchy- Schwartz, putting g in one term, / in the other term, and half of 
all the rest in both terms. Thus ((26|) follows with the elementary inequality that 
we just used in the previous estimate. 

To obtain ([?5|) . on the other hand, we will bound above the r.h.s. of the last 
displayed inequality by splitting into the regions jw — < 1 and |t; — ti^j > 1. On 
the former region mI^^ < , Vm > so we may take the L°° 

norm of (v) ™ use (|27[) . and apply Cauchy-Schwartz as above to obtain 

<2'"'\h\^.^Jg\L2^Jf\Li^, Vm>0. 

On the latter region, |u — u*| > 1, we have Ml^'^\v — v.^,p^'^'^ < Af^^ {vy^'^'^ . 

Applying Cauchy-Schwartz as above in both the v and integrals gives 

Summing these two inequalities gives psp . □ 
Proposition 2. For all ^ e R, m > and integers k, the inequalities are uniform: 
(29) |T^^(g,/.,/)| <22^'=|5|^k.|^«^/^|^^,J^V|l^^,^, 



(30) K\9.h,f)\<2'-^%\L.__Jw'h\^.>^ Iz^VIl 



1,1 + 23 



Moreover, 

(31) |ry(5,0,/)| + \Ty{gJ,^)\ < 2'^%\L2_Jf\, 

where as usual (p is defined as in (j23[) . 



Proof. As in the previous Proposition [l] the key to these inequalities is the sym- 
metry between h and / coupled with two applications of Cauchy-Schwartz. The 
difference is that, this time, the Carleman representation ((2T|) will be used and the 
main integrals will be over and v' . In this case, the quantity of interest is 

|2 _ I,, _ ,/ 



b [ !L !*' , \l ) r rL-9 Xki\v - v'\). 



The support condition yields \v ~ v'\ ^ 2 ^. Moreover, since &(cos0) vanishes for 
9 G [7r/2,7r], we have \v' — u*| > |w' — v\. Consequently, the condition (jSj gives 



\v' — v.f\^ -\- \v — V \^ I ■ \\v' — vA"^ 



Thus, the integral is bounded by a uniform constant times 



|n-l+2s 



I El' I"" - V' 

As a result 



T 7iLi+2. 1^^ - y*r'^hk{\v - v'\) < 2^-^'^\v' - 



Ty{g,h,f)\ <22^^- / dv' / dv,, 014 \v' ~v,p+^'\gMf \ w^\v'). 
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Now the relevant estimates can be established as in Proposition [T] 



□ 



We will now estimate the operator T^'^, which is more difBcult and more technical 
because it contains the post-collisional velocities ([2]). A key problem here is to be 
able to distribute negative decaying velocity weights among the functions g, h, 
and /. In the previous propositions this distribution could be accomplished more 
easily. Because the weight is in the v variable and the unknown functions are of 
the variables v' and v'^^ we have to work harder to distribute the negative weights. 
These decaying weights are used to close our estimates for the soft potentials. In the 
following we generalize methods from |35[l36j to overcome this particular difficulty. 

Proposition 3. Fix an integer k and £'^,i~ > 0, with i = £^ — £~ . For any 

< i' < , with i + £' ^ £^ — {£^ — £'), we have the uniform estimates 



(32) 



\w'+''h\LAw'f\L^^ 



We also have a similar estimate with the roles of g and h reversed 



(33) 



Tl'Hg.hJ) <2^^^\w'^-''g\LAh\^K>^ Wf\ 



i + i',-, + 2 

e+e' 



Froof. We plan to estimate T^' (p, /i, /) in ((20)) from above as 



< 



dv 



dvt 



da Bk{v - v,,a) w'\v) \g',h'\ \f\ ^ M.. 



We will utilize a number of splittings of this upper bound. Firstly, when > ^\v\ 



we have w^^{v)mI''^ < [M^MY for some 5 > 0. Then T_f:'^(.g, h, f) 
is bounded above by a uniform constant times 



on this region 



dv 



da B,iv-v,,a) \g:\\h'\\f\ M^M' 



Then, with Cauchy-Schwartz, this quantity is again bounded by a uniform constant 
times 



dv / dv.^ 



da Bk{v-v,,a) \g'J'\h'\' {MiM')'^' 



1/2 



dv 



dv, daBk{v-v,,a) {M,Mf^^\f\^^ 



1/2 



For the second term, with (|28p . we clearly have for any m > 



dv 



dv, / da Bk{v - v,,a) {M^uf^ \f\^ < 2^'''\f\i 



For the first term, do another prc-post and take the L°° norm of M^^^g, to obtain 



that 



T, ■ (g, h, f) , when restricted to this region, with ([27]) is bounded by 



-)2sfc 
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(times some uniform constant for any fixed m). If we instead took the L°° norm of 
M^/'^h and applied the same procedure, we would have the same inequality except 
with the roles of g and h reversed. 

Our next splitting is onto the region \v^,\ < ^\v\. Notice that if |w — w»| < | 
then > — |w — w^j > ^\v\. Thus |w — w*| > ^\v\ when < ^\v\. Further 
suppose that \v\ > 1. On this region, from ([2]), we have, in particular, 

\v -v.\> ^\v\ > (1 + + \v\'y/' > (1 + + 



T!:'\g,h,f) 



re- 



This lower bound will be used several times just below. Then 
stricted to this domain is bounded with Cauchy-Schwartz by a uniform constant 
times 



dv I dv^ / da Bk{v -v^,a) w^'iv) |glr|/iTvA/, 



dv I dv, da Bkiv-v,,a) \f\^ w^\v)^/m. 



1/2 



1/2 



(34) 



The estimate for the second term involving / follows as in (|28p . The estimate for 
the first term involving g'^h' requires more explanation. First of all, on the region 
of interest using also the collisional geometry ([2]), = w^^*" which is 

2e/ \ ^ 21'^ I i\ 21+/ i\ti I I 'I 1 I ^ 2U+-£')/ r\ 2(e+e')/ i\ 

W (V)^VU (V )w (^'*) (1 + P I + ^ (V,)W ^ (^ )• 

With that estimate, we take a pre-post collisional change of variables, estimate 
< 1 and \v — v^p^'^'^ < {v^y'^'^^ when 7 + 2s < on the region of interest. This 
combined with ([28l) yields the estimate above. 

Finally when \v\ < 1 and \v^,\ < ^\v\ we have that jw^p + \v\'^ < 2. This region is 
invariant under the pre-post collisional change of variables. We can thus generate 
(artificial) Maxwellian type decay in both variables and v since all the variables 
are trivially bounded. In this last region, with these preparations, the estimates 
follow exactly as in the first splitting which was \v,\ > ^\v\. Thus when 7 + 2s < 0, 
we obtain stronger results than ([5^ and ([33]) (and may, in fact, omit the second 
line of terms on the right-hand side for both inequalities). 

In the case 7 -I- 2s > 0, we use a particular sequence of Cauchy-Schwartz inequal- 
ities which has already found utility in the cut-off context by Guo [37] : 

..21/ 



dv{ I dv, I daBkiv^v^a) \g',h'\ ^ M, j w''{v) \f\ 

1/2 / » N 1/2 



< Jdv(^J dv,da Bly^}j (^j dv,da \g',h'f^/M^^ uj^'^{v) \f\ 

1/2 

< 2^^'= J dv {vy+^' (^j dv,da \g',h'\^ w^'{v) \f\ 

<2'^'^\w'f\,2^^^ dvdv,da {vy^''w'^''~'\v:)w'^'+'Hv')\g:h'\'^' 
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We have also used w'^^{v)\/M^ < w'^^^^ ^ {v'^,)w'^^^'^^ H''^')^ which follows from the 
splitting as in the previous case. Since {vy^'^'^ < {v'y^'^" + (w*)'*'^^* , we have 
established the estimates ((32|) and ([33]) for 7 + 2s > (the additional terms on the 
right-hand sides of ([32]) and ((33|) arising from the inequality just observed). □ 

Proposition 4. Suppose that k > 0. Also consider £^,£~ > 0, with I = — t~ 
and < i' < £^ so that i -\- 1' = ~ {(T — £'). We have the uniform estimates 



(35) n'^g^hj) < 2'^'^\g\^K, h\^.^Jw'f\^ 

We also have 



(36) 



Tl^\g,hJ) <2'^'\w'^-''g\L.\h\ k'/lL^, 



Proof. Notice that both inequalities are proved exactly in the proof of Proposition 
[3] in the case when 7 + 2s < 0. For 7 + 2s > 0, the only estimate which requires 
modification in the proof of Proposition |3] is the one in the splitting surrounding 
the estimate (p4|) . This estimate can be modified as follows. When k > 0, we 
have |w» — < 1 which implies that for some S > there is exponential decay as 
< {M^M'y. This decay further implies 

Now in the estimate for (|34|. when 7 + 2s > 0, this polynomial growth can be put 
entirely on h' , which finishes the argument. □ 

This concludes our size and support estimates. In the next sub-section we will 
prove estimates which incorporate the essential cancellation properties of the Boltz- 
mann collision operator in the appropriate geometric framework. 

2.3. Cancellations. We now establish estimates for the differences y^'^_ j"^'^ 
rpk.e _rpk,e ([20]) and ((2T|) . We wish the estimates to have good dependence on k 
(the estimate should contain a negative power of 2'^), but this improved norm will be 
paid for by assuming differentiability of one of the functions hov f . The key obstacle 
to overcome in making these estimates is that the magnitude of the gradients of h 
and / must be measured in some non-isotropic way; this is a fundamental point, 
as the scahng is imposed upon us by the structure of the "norm piece" {Nf, /). 

The scaling dictated by the problem is that of the paraboloid: namely, that the 
function f{v) should be thought of as the restriction of some "lifted" function F 
of n -I- 1 variables to the paraboloid (?;, i|wp). Consequently, the correct metric 
to use in measuring the length of vectors in R" will be the metric on the parab- 
oloid in R""""^ induced by the (n -I- l)-dimensional Euclidean metric. To simplify 
calculations, we work directly with the function F rather than / and take its (rt-|-l)- 
dimensional derivatives in the usual Euclidean metric. This will be sufficient since 
our Littlewood-Paley-type decomposition in Section |3] will grant a natural way to 
extend the projections Qjf into n + 1 dimensions while preserving the relevant 
differentiability properties of the n-dimensional restriction to the paraboloid. 

To begin, it is necessary to find a suitable formula relating differences of F at 
nearby points on the paraboloid to the various derivatives of as a function of 
n + 1 variables. To this end, fix any two v,v' S R", and consider 7 : [0, 1] — t- M" 
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and 7 : [0, 1] 



pri+l 



given by 



7(61) = 6v' + (1 - 6)v, and 7(6*) = [ov' + (1 - e)v, ^ \6v' + (1 - e)vf 

Note that 7 hes in the paraboloid {(wi, . . . , Wn+i) £ M"^"'^ | Un+i = \Wi + ■ ■ • w^) }, 
and that 7(0) = u and 7(1) = 1;'. Also consider the starred analogs defined by 

7*(^) = + and 7*(^) - + (1-0)1'*, i|0< + (l- 0)i;*|') . 

Elementary calculations show that ^{9) + 7*(^) — v and 



^ = (t;'~i,,(7(0),z;'-i;)) and ^ = (0,1^;'-^'), 



^7* 



— (u' — w, (7,(6*), w' — u)) and 



= (o,K-H')- 



Now we use the standard trick of writing the difference of F at two different points 
in terms of an integral of a derivative (in this case the integral is along the path 7) : 



F v\- 



n'|2 



- F V 







(37) 







' de(^^-{\7F){2{e))], 



where the dot product on the right-hand side is the usual Euclidean inner-product 
on M"+^ and V is the (n + l)-dimensional gradient of F. For convenience we define 



\\/\'F[vi,. . . ,Vn+i) = max sup 

0<J<» |5|<i 



(i-Vy F{vi,...,Vn+l) 



1 = 0,1,2, 



where ^ G and |^| is the usual Euclidean length. In particular, note that we 

have defined |V|°i^ = 

If V and v' are related by the collision geometry ([2]), then {v ~ v' ,v' ~ v<,) = 0, 
which yields that 

{-f{e),v' -v) = {v,,v' ^v)^{i^e)\v~v'\^, 

(7*(6l), u' - w) = (f*, w' -v)+ 9\v - v'\^. 
Thus, whenever — w'| < 1. which holds near the singularity, we have 



<\v - v'\ (w*) . 



Throughout this section we suppose that |w — u'| < 1 since this is the situation 
where our cancellation inequalities will be used. In particular, we have the following 
inequality for differences related by the collisional geometry: 









dO 


+ 


dd 



(38) 



F v' 



„'|2 



F V 



< {v,)\v^v'\ / d9 |V|F(7(0)). 



Furthermore, by subtracting the linear term from both sides of ([57)) and using the 
integration trick iteratively on the integrand of the integral on the right-hand side 
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of ([37|). we obtain 

l,/|2 



(39) 



F v' 



-F[v,^]~%0)-^F{v) 



dO 



< 



{v.Y\v-vr deM'F{^{e)). 



We note that, by symmetry, the same result holds when the roles of v and v' are 
reversed (which only changes the curve 7 by reversing the parametrization: 7(0) 
becomes 7(1 — 0)). It is also trivially true that the corresponding starred version 
of (|39p holds as well. We will use these two basic cancellation inequalities to prove 
the cancellation estimates for the trilinear form in the following propositions. 

Proposition 5. Suppose f is a Schwartz function on M" given by the restriction 
of some Schwartz function F on R"+-'^ to the paraboloid {v, Let |Vp/ be the 

restriction of |Vpi^ to the same paraboloid. Then, for any k > 0, 



(40) 
(41) 



{t1'' - T'i''){g, K f) < 2^''-'^'\g\^Kf^ \w'h\ 



wWf 



Each of these inequalities hold for any m > and any £ CzR. 
Proof. We write out the relevant difference into several terms 

'dj, 



dj* 



dO 



(0)- VAf^(w,) 



+ {w^\v') f - w^\v) f) -=(0) • VMpM 



d9 

+ MpM (^w^'iv') f - w^'iv) / - J(0) • ^iw^'f){v) 

+ Mpiv,) (^(0) • v{w^'^f){v) \ = I + n + in + IV + V. 



We further split (r_^'^ - T'l'^){g,hJ) = + T° + T™ + r"' 
corresponds to the first term in the splitting above, etc. 
We begin by considering the first and last terms, that is 



T^, where 



rpH 

and 



dv / dv^ 



^7* 



da Bk{v - «,,a) g,h \ -={Q) ■ ^Mp{v,) ) w''{v) f, 



1 dvl 


dv, [ 




i" Js 



da B,iv~v,,a) g,h Mp{v,) ( ^[Q) . ^ (w'' f){v) ) . 



^7 



We may estimate both of these terms in exactly the same way. 
First we define the extension. If / extends to R""*"^, then 



{U?'^ f){vi, . . . ,Vn+l) = W^^(w)/(wi,...,U„+l) 



21 1 
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Here we think of w^^(u) as being constant in the {n + l)-st coordinate With this 
extension it follows that for i = 0, 1, 2 we have 



(42) 



21 \ 



We will use this basic estimate several times below. 

For , notice that ^(0) is linear va. v' ~~ v and has no other dependence on v' . 
Thus the symmetry of Bk with respect to a around the direction |^~^* forces all 
components oi v' — v to vanish except the component in the symmetry direction. 

\v-v'\'^ 



Thus, one may replace v^v' with "_^^*| 



in the expression for -^(0) 



Since (v 



V , V 



0, the vector further reduces to 

,'|2 



\v — V 



The last coordinate direction of 



\v — 
US 

l\2 



< 2 



-2k \ 



-1 



Hence 



(0) is given by {v, v' — v) which reduces to 

< (l + \v-vJ~') 



These bounds allow us to employ the same methods as used previously. To be 
precise, we must control the following integral 



-2k 



dv 



dv^ I da Bk Afi 
/s2 



Here we absorb any powers of (w*) by for any small e > 0. The estimates pOj) 
and (PT|) for then follow as in the proofs of and (psj) . the only difference 
being the extra factor jw— |~^; however, since 7+2s— 1 > — (71— 1) one may replace 
7 with 7 — 1 in the earlier proofs without fear of destroying the local integrability 
of any of the singularities arising therein. 

The estimation of can be handled in exactly the same way because -^(0) also 
depends linearly on v — v' and has no other v' dependence. Here we have exactly 
the same estimates for -3^(0) as just previously obtained for -^{0) by symmetry. 
Then (gOl) and jH]) for follow again as in the proof of ^ and ([25]). As in (|42|) . 
we use 

\V\'Mf3<y/M, i = 0,1,2, 

where the extension is defined as Mp{v) = (27r)~"/^ P/iivi, ■ ■ ■ ,Vn) 6"'""+^/^. 
We now turn to the estimation of the term T™, which can be written as 



T™= / dv dv, da Bk{v-v,,a) g,hw^^{v') f 

/ , ^7* 
X (MM) - MpM - ^(0) • VMp{v,) 

We apply the starred analog of (p9| to obtain 

da Bk\g.h\w'''{v')\f'\{AUACr, 





dv 1 


dv, / 




I" Jl 


i" Js 



where we have just used the estimate ^ M{'^,{9)) < {M,MlY, (valid for all suf- 
ficiently small e > 0) which follows directly from (v,) < {j*{0)) < (w*) since 
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w* = 7* (6*) + 0(v* — v'^) and |u* — < 1. At this point, after applying Cauchy- 
Schwartz, it suffices to prove two estimates. Tlie first is that for any 7ti > 



(43) 



dv / dv, / daBkMl w^^{v) \g^'^\h\ 



\g\L^Wh\i 



where ur {v) replaces ur {v') because (w) « {v') when \v — v'\ < 1. The second is 



dv 



dv^ 



daBk{M,M',Yw^\v')\r 



pl\2 



< 2" 



'f\i 



Both of these inequalities follow exactly as in the proof of and 

It remains to only prove (pOj) and (|^T|) for T"^ and T"'^, both of which involve 
differences in w'^^f. We use the difference estimates (|55)) for T^" and ([55]) for 
to obtain that for any fixed e > both |r°| and |T^^| are controlled by 

<2-2'= / ^6* / / dv, I da Bk Ml'^ w'^\v) \g,h\ |Vp/(7(f))- 



The loss of e comes from the factor (w*) in ((38|) and ([39]) in addition to the bound 
|V|M^(-(;*) < mI and the estimate for ^= . These also account for the 2^^*^. 

Note, though, that the factor 2~^'^ comes directly from (p9| . but (|38|) only furnishes 
a factor of . In this case, there is an additional factor of 2^'^ available in the 
estimate for T^" arising exactly from the derivative of 7*. Also, {v) « (which 
accounts for the replacement of w'^^{j{d)) by w'^^{v)). 

With that last estimate and Cauchy-Schwartz, it suffices to use (^5)) and show 



(44) 



d9 



dv / dv. 



daBkMl-'w^'hie)) |V|V(7(^)) 



<2^'=|w^|Vp/l 



This uniform bound follows from the change of variables u — 9v' + (1 — 9)v, which 
is a transformation from v to u. In view of the collisional variables ([2]), we see (with 
5ij the usual Kronecker delta) that 



dui 



dv' 



dv^ 

with the unit vector k - 

duj 



dvi 



dvi 



(u — w,)/|u — |. Thus the Jacobian is 



(k.a) 



Since h{{k,a)) = when (fc,cr) < from and G [0,1], it follows that the 
Jacobian is bounded from below on the support of the integral (|44| . But after 
this change of variable the old pole k = (w — )/!?; — f*| moves with the angle 
cr. However when one takes k = {u — w*)/|u — then 1 — (fc, cr) « 1 — {k,a) , 
meaning that the angle to the pole is comparable to the angle to k (which does not 
vary with a). Thus the estimate analogous to ([^5)) will continue to hold after the 
change of variables, giving precisely the estimate in (j^J). □ 
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This completes our proof of the cancellations for the a representation as in ([20| . 
In the following, we estimate the cancellations on h instead of putting them on /, 
and for this we use the Carleman representation as in (j21[) . 

Proposition 6. As in the prior proposition, suppose h is a Schwartz function on 
M" which is given by the restriction of some Schwartz function in W'^^ to the 
paraboloid {v, jl^P) and define |V|'ft. analogously for i = 0, 1, 2. We have 



(45) 
(46) 



\w'f\i 



k,e\ 



The above inequalities hold uniformly in k > 0, m > 0, and £ G 



Proof. This proof follows the pattern that is now established. The new feature in 
(|45)) and (|46l) is that, from ((2T|) . the pointwise differences to examine are 



Mp{v[)h^{v - v,)\v -v.r- Mfi{v,)h'<^{v' - v,)\v' - v,\' 



de 



(0) • VM^(w,) h'^{v' - v,)\v' - 



+ (^^(0) ■ VM/5(w,)j ih<^>{v - v,)\v -v,\"~ h'^v' - v,)\v' - 

+ (^Mpiv'J - Mp{v,) - ^(0) • VM^(w,)^ h<^{v - v,)\v - 

+ Mp{v,) [h^{v - v,)\v - «,|" - h'^{v' - v,)\v' - - j(l) • VF(7(1)) 

+ M;3(«.)(^(1)-VF(7(1)) 
Here since jg\l{()) — w»P = 2 {v' — v,j{9) — v^) = 0, the canccUation term is 
J(l) • VF(7(1)) = <i>{v' - v.)\v' - v.r^il) ■ Vh'. 

Wc again split (T_f;'^-T^^)(g, h, f) = Tl+Tf+T™+Tl^ +T"^ , where corresponds 
to the first term in the splitting above, etc. 
For the last term Tj, we have 



rpV 



dv' [v') I dv^ / d-Kv — 
Jk" Jei:' ^ 



Bk 



I?;' — t!*!!?; — w*!" ^ $(u — w*) 



de^ ' 

In this integral as v varies on circles of constant distance to u', the entire integrand 
is constant except for -^{X)- If we write -^{X) ^ sum of two vectors, one lying 
in the span of the first n directions and the second pointing in the last direction, it 
follows that we may replace the former vector by its projection onto the direction 
determined by f' — w*. But since the original vector points in the direction v — v', 
the projection vanishes. Since the last direction of ^(1) is exactly {v' , v' — v), the 
corresponding integral of this over v also vanishes by symmetry. 
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Similarly, we consider the first term T^, we have 



dv w [v ) I dv^, / dTTy— — — ; — -—, r 



Here the explanation is the same as in the previous case. 

The remaining terms incorporate cancellations. To estimate T}^ and Tj^ we 
study the function of 9 which measures differences given by 

h{2{e)Mm-v,M9)-v,r. 

Notice that both terms II and W include the difference in the values of this function 
at 6* = and 9 = 1 (we also extend h to R"+^ and make the de facto extension of 
$(w — — to R"+^ by taking the new function to be constant in the last 
variable) . In terms of II and IV, our operator from ((2T|) takes the form 

^11 ^^iv ( ^ I 21, [ , f , Bk gj' (II +IV) 

We use ([38|l and p9|) to estimate II and IV. The additional estimate needed is that 

V''($(« - v,)\v - v,n\ <{l + \v- v^Di^v - v^)\v - 

which simply comes from differentiating with respect to v (since the extension is 
taken to be constant in the last direction, the gradient V reduces to the usual n- 
dimensional gradient). Applying ([38)) (with the estimates for ^ ) and (p9)) gives 
the estimate 

Again, |ti — < 2~'^ and also |w — w |7(^) — ^'*|- With all of that we may 
estimate the terms iTjJ^I and ItJ^I above by the following single term: 



-2k 



'd9f dv' f dvj dn.. . . B,{v- v^,2v'~v^v.) 



JS." JR" J E. 



2 min{|i; — 1} 
X w^'W) mI%J'\ |Vp/i(7(0)). 



The estimates required now for the term above are completely analogous to a cor- 
responding estimate from a previous proposition. In particular, we split into the 
regions |f — < 1 and |w — > 1. When |u — > 1 we apply Cauchy-Schwartz 
to the quantities above carrying \g*f'\ into one of the integrals and \V\'^h{'-^{9)) 
into the other integral. The rest of the terms above are carried half into each in- 
tegral. The integrals involving |Vp/i(7(6')) arc estimated as in (|44|) after changing 
back from the Carleman representation to the a representation; see for instance 
[551 Proposition 10]. The integrals involving \g*f'\ are estimated as in the proof of 
([25[) . ([25[) . and Proposition [1 

When instead |u — w.^! < 1, since also |u — w'| < 1, we have strong exponential 
decay in all variables. After taking the L°° norm of either (7(6*)) ~™ |Vp/i(7(6')) or 
(z;*)""* and using ([77]) . we can use the Schur test for integral operators to finish 
off our estimates, or simply Cauchy-Schwartz as in Proposition [1] and Proposition 
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[21 The additional singularity does not cause a problem because 7 + 2s — 2 > —n 
which implies local integrability. 

Regarding the estimate for T™, we notice that the estimate of holds if v, v' 
are replaced by w*, v'^ and we replace 7 with 7*. Using (j39]) in this case, we have 

|Tf I < 2-^^ de I dv' I dv. 1^ ^ dTT,, Bu{v~v., 2v'~v -v^.) 



X A/y2|5,/'| \h\. 



Now this estimate can be handled as in the previous case for T^^. This case is easier 
because there is no 7(0) here, which required another change of variable. □ 

This concludes our cancellation estimates for the differences involving three ar- 
bitrary smooth functions. We will also need cancellation estimates when have a 
more specific smooth function satisfying the following estimate 

(47) |Vp0< C^e-=l''l', C^>0, 00. 

Above is any Schwartz function on M" which is given by the restriction of some 
Schwartz function in M"+^ to the paraboloid (w, ^ji"!^) and |Vp0 is defined analo- 
gously as usual. With this in mind, we have the next estimates: 

Proposition 7. As in the prior propositions, with (|47p . for any k > we have 

{Tl''^Tl-')ig,h,4>)\ + \{T'l'' -Ty){g,cb,h)\<C^ 2(^''^-''>%\L._Jhy^. 

The above inequalities hold uniformly for any m > and £ G M. 

Notice that the proof of this Proposition [7] follows exactly the proofs of the 
previous two Propositions [5] and [HI The main difference now is that when going 
through the last two proofs above, as a result of (|47l). we will always in every 
estimate have strong exponential decay in both variables v and w*. This strong 
decay allows us to easily obtain Proposition[7]using exactly the techniques developed 
in this section and the last one. We omit repeating these details again. 

2.4. Compact estimates. In this sub-section we prove several useful estimates 
for the "compact part" of the linearized collision operator ([T3l) . Here we use the 
integer index j instead of k to contrast with the kernel k below. Our first step is 
to notice that the Carleman representation ([?!]) of Ti^^{g, Af^^, /) grants 

T^\g,Mf,,J)= J^^dv' w^\v') f{v') J^^dv, gM k]+'\v',v.), 
where for any multi-indices /3 and /3i we have the kernel 

Recall the domain = {v' E R" : {v^ — v' ,v — v') = 0} , and that diTy denotes 
the Lebesgue measure on this hyperplane. More generally, suppose that 

where satisfies |0(w)| < C^e"^'"' , and similarly < C^e"^''"' for any posi- 

tive constants C^, C^, and c as in ([23]). Then we have the main compact estimate: 
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Lemma 7. For any £ G M, Kp'^ satisfies the following uniform in j < estimate 



dv' 



Furthermore, the same estimate holds if the variables are reversed 



dv^, 



The crucial gain of the weight — (rt — 1) has been known in the cut-ofF regime; 
see Grad [32] , and also [15] . The essential new feature here is to gain this weighted 
estimate without angular cut-off using our decomposition of the singularity and 
various estimates involving ([49| below. 



Proof of Lemma [^ The second upper bound follows easily from the first after notic- 
ing that the roles of v' and can generally be reversed in the first upper bound 



for 



given below in 



We focus here on the first bound from above. We perform the change of variables 
= V — v' in V, with as usual v' = v + — v' , then the kernel can be written 



d-K 



Bj {v' — + z,v' ~ — z) 



+ v^)(f>(z + v'). 



Above uj = |"^,_^* and u}^ denotes integration over the n— 1 dimensional hyperplane 
orthogonal to w with Lebesgue measure dir^ ■ We expand the exponents as 



\z+v^\^ + \z+v'\^ = \z+v^\^+\z+v^ + {v'-v^)\^ = ]^\v' -v^+2{v^+z)\^ + \\v'-v^\^ 



Furthermore, the kernel can be expanded as 



Bj{v' — v^ + z,v' —v^: — z) = b 



w' — w* + z 



V — V.J, — z 



+ zY \v' — — z| 



<^{\v'-v,+z\)xM- 



The orthogonality condition implies ( {^J^Xl^ {irZ-l] ) = ■ Then, 



since b{t) < {I - ty 



from ([3]), we have 



w — W» + Z V — v^, — z 



\v' — u, -I- z| ' \v' — — z\ 
We thus conclude that 



< 



\v' — -I- zp 



(48) nf^{v',v.) 



< 



exp 



diTz \v' — + z| 



y+2s+l XjXkl) -2c\z+ 



|n-l-|-2s 



We will estimate this upper bound. Without loss of generality, in the following for 
notational simplicity we set the unimportant constant to c — 1/4. 
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Suppose 7 + 2s > — (ji — 2). For the integral over z e w^, with j < 0, we have 

dn, \v' - + z|T+2^+i 



|^|n-l+2s ^ 



|n-H-2s 



|2|ri-l+2s 

We have used that |z| > 1 when j < 0. The upper bound just above is 

dn, {v' ~v, + z-{v,-{v,,Lu)Lj)y+^'+' e-l^"l'/2 

^ Xj(l^- (w* - (w*.^)^) I) 



\z- (u* - (u*,a;) |"-i+2s ■ 
Now we have the uniform, pointwisc inequahty (vahd for any p e M) 

With that, we observe that the previous displayed quantity is bounded above by a 
uniform constant times 

/ / \ \7+2s+l / / i|7+2s+l| 



X / dn^ 



Xji\z - jv^ - {V^,UJ)UJ) I) _|j|2/4 

The essential new claim is the following uniform estimate: 



(49) 



\z~ {V^ - {v^,Uj)uj) |"-l+2.s'^ 



Modulo (|49| ; the combination of the estimates above shows that 



\V' — V, 



exp 



^ , / \ \7+2s-(n-2) / / i|7+2s+l 



We thus deduce that 



Kf^{v',V,) W'{v')dv' 



dv' 

< C^C^ 2"'^ I dv' {v, - {v,,uj)ujy+^'~^'''^^ \v'-v,\-^w'{v') {v' - 



X exp 



This leads directly to 
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with 



J' = 



/ \ \7+2s-(n-2) / / \|7 



+2s+l\ - W*) 



X exp 



\v' ~ f*| 



Notice that up to a change in parameters this is the same J' which appeared in 
Mouhot-Strain [49l p. 525]. Therein the fohowing basic pointwise estimate is shown: 
J' ^ (w*)''^^''^^""^^ This completes the estimate, once we prove p9)) . 

We now estabhsh the key estimate (|49| . Since j < and Xj is supported on 
p-i-i^ 2"-'], we see that \z — (u* — (w*,a;) w) | « 2"-' > 1. Thus in particular, since 
also |xil-L°° ^ Ij 'we have the following estimate 



|z- (u* - (u,,a;) |«-i+2s 



/4 



This concludes our proof of the new claim in ()49|) . and the proof. 



□ 

With the compact estimate for the kernel from Lemma [7] in hand, we can auto- 
matically prove the main estimate for the compact term (|50p below. Using similar 
methods, but without Lemma [71 we will also prove ((5T|) and ([52]) . 



Proposition 8. Consider (f> satisfying (|23p . FFe have the uniform estimates 



(50) 



< 2 



2sfc L,,« 



These hold for any k < 0, £ £ W. Additionally for any m > and any k we obtain 



(51) 

Furthermore 
(52) 



T^''(5,/,0) 



2 



2sfc 



\9\L^Jf\l 



Tl:'\<t>,hj) 



i?ac/i of these estimates hold for all parameters in ^ and 

Proof. For ([50]) . using the formula in ([2T|) and Cauchy-Schwartz we have directly 

dv^ \g{v<f)\ Kf''^{v',v^) 

1/2 



< 



< 



^d^;' 1.2^(1;') |/(«')| 
dv' w'\v') \f{v')\'' 

dv^ |5(w*)|^ 



dv 



Kf''''{v',V^) 

dv' w^\v') 



iv',v*) 



1/2 



Above we consider the kernel k^''^{v\v^) with t/i = Mp and as in Now we 

observe that Lemma [7] immediately implies (|50p . 

To prove ([5T|) . notice that the bound ([23]) implies there is strong exponential 
decay in both variables v and f*. The estimate (|51|) then follows from the Schur 
test for integral operators or simply using Cauchy-Schwartz as above. 
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The proof of ((52|) follows from a similar application of Cauchy-Schwartz. The 
main difference is that this time we have strong exponential decay in both v'^ and 

which displays the appropriate symmetry under the needed prc-post collisional 
change of variable. □ 

This completes our basic compact estimates. In the next section, we develop the 
geometric Littlewood-Paley theory adapted to the paraboloid in n + 1 dimensions. 



3. The ?7.-dimensional non-isotropic Littlewood-Paley decomposition 

In this section, we expand the development of the non-isotropic Littlewood-Paley 
projections that we initiated in [33]. The principal new result that is needed in this 
direction is an understanding of the interaction (specifically, the commutator) be- 
tween the geometric Littlewood-Paley projections and the isotropic velocity deriva- 
tives dp . These additional estimates are necessary because of the higher derivatives 
present in the norm N'^'^, which are ultimately necessary because the singularity 
of the kinetic factor is strong enough that many of the estimates in [321 must 
be replaced with L°° estimates. These L°° estimates, in turn, are related back to 

based spaces via Sobolev embedding. 

Throughout this section, wc will use the variables v and v' to refer to independent 
points in M", meaning that wc will not assume in this section that they are related 
by the collisional geometry. The reason we choose to use these variable names is 
that they give a hint about where the Littlewood-Paley projections will be later 
applied in situations which do involve the collisional geometry explicitly. 

3.1. Definitions and comparison to the non-isotropic norm. Following the 
procedures and notation of [33], we first describe the lifting of vectors v e M" to 
the paraboloid in W'+\ Specifically, for v € R", let v = {v, ^\v\^) e and 
consider the mappings : R" — )• R" and t^, : R" — )■ R"+-'^ given by 

TyU u - {1 — (v) ^) {v,u) \v\^'^v, and Lv'^'= ('''v^A''^) ^(^'j""))- 

These mappings should be thought of as sending the hyperplane = to the 
hypcrplanc tangent to the paraboloid {v, ^|wp) at the point v- It's routine to check 
that (v^Tyu) = (v) ^ {v,u) and |T„up ~ |up — (v) ^ (u, li)^. which implies 

meaning that is an isometry from one hypcrplanc to the other. Moreover, it is 
not a hard calculation to check that 

1, ,2 

V + TyU = V + TyU + - \u\ e„+i, 

where e„+i == (0, ... ,0, 1). Next, fix any C°° function supported on the unit ball 
of R"+^ and consider the generalized Littlewood-Paley projections given by 

(53) P.fiv) f dv'2-=^{V{v - i/)) {v') f{v') 
Q-jf{v) = P,f{v)-P,^,f{v), J>1 

(54) = ( dw'2"^V'(2^'(w-2L)) 
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where ip{w) = ip{w) — 2~^ip{w/2). Here Pj corresponds to the usual projection onto 
frequencies at most 2-' and Qj corresponds to the usual projection onto frequencies 
comparable to 2^ (recall that the frequency 2^ corresponds to the scale 2~^ in 
physical space). We also define Qo "= Po to simplify notation. In order for these 
projections to be generally useful, the function ip must be chosen to satisfy various 
cancellation conditions which we will discuss later. For now, we note that, as long 
as the integral of Lp over any n-dimensional hyperplane through the origin equals 1 
(which will be the case for any suitably normalized radial function), we have that 
-Pjfi'v) ~^ fi"^) &s j ^ oo for all sufHciently smooth / and that 



(55) 



dv\p,f{vr{vy 



< 



dv\fiv)f (v) 



uniformly in j > for any fixed p G K and any p £ [1, oo) (with suitable variants 
of this inequality also holding for p = oo as well as for the operators Qj ) . 

The principal reason for defining our Littlewood-Palcy projections in this way is 
that the particular choice of paraboloid geometry allows us to control the associated 
square functions by our non-isotropic norm. This informal idea is made precise in 
the following proposition (the proof of which may be found in [33]). 

Proposition 9. Suppose that \Qj{l)(v)\ < 2~^^ holds uniformly for all v € M" and 
all j > 0. Then for any s G (0, 1) and any real p, the following inequality holds: 



(56) 



3=0 



2s] 



dv \Q,f{v)\''{v 



I/I 



LI 



+ dv dv {{v){v)) ld(^^y 



d{v, v') 



')<!• 



This is true uniformly for all smooth f . 



3.2. Littlewood-Paley commutator estimates. For convenience, let us abbre- 
viate 2"-'(/3(2^w) fjiw) and likewise for ipj. We seek at this point to relate the 
corresponding Littlcwood-Paley square function to the norm Nf '"^ in arbitrary di- 
mensions 71 > 2. This was done in three dimensions for N"''^ in our first paper [33) . 
For this extension, it will also be necessary to make estimates involving the commu- 
tators of the Littlcwood-Palcy projections with isotropic, n-dimcnsional derivatives 
(i.e., derivatives which do not respect the intrinsic geometry). The commutators 



themselves. 



Qj fiv) = [Qjf] (v) - Qj -irj (^)' are easy to calculate 



d_ 

dvi 



,Q3 



f{v) = I dv'2' (V.l^ + V,Vn+lA {v - z/) {v') f{v') 

.Ion \ / 7 
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Above Vi is the i-th component of V. After an integration by parts, 
d 



dvi 



, Qj 



f{v) = (V,i^ + i..V„+iV')^. {v - 2/) {v') f{v') 

dv'-^ [<P,{v^£) {v')]f{v') 



dv'2'{v, - v[) (V„+i^j^ [v-^) {v') f{v') 
dv'tpj{v^2/_)vl{v'y^ f{v'). 



In particular, this commutator may be written as 

d 



dvi 



f{v) = Q,f{v) + Q,f{v), 



where Qj is given by replacing 'ip{w) with WiV n+ii^iw) in the integral (j54p and 
f{v') == v[{v')~'^ f(v'). The end result is that, after taking /3 derivatives and 
studying these commutators, we may always write c?^2~l"l^ V"Qj as a finite sum 

l/3l|<l^l k 



where each has a form the same as (|54)) for some V' , the derivatives (3i satisfy 
^ and the weights cjfc(u') are polynomials times powers of {v') which 
together tend to zero at infinity. Above V represents the n + 1 dimensional gradient 
which is restricted to the paraboloid, (u, P), post differentiation; since this is 
the way in which that particular operation is used in the sequel. 

We are able to compare these weighted non-isotropic Littlewood-Paley projec- 
tions to the non-isotropic norm using the same method which we used to prove 
Proposition El in [33], that is, by completing the square. We bound above the 
expression 

dv [ dv' I dzi;''Av-z)^){y^-i){f{v)- f{v')f {zY {v){v')LOk{v)u:k{v'), 



by integrating over z and comparing this to the semi- norm piece of N^'^, then we 
show that this term is also equal to 



dv |QfK/)(z;)|2(«) 



dv Q){uJk){v)Q){u:kf){v) {v)" , 



(in both expressions, we suppressed the dp-^^ acting on /). As long as one has the 
following uniform inequality for all j > and all u, 

(57) \Q^iu;,)iv)\<2-'^, 

these estimates may be multiplied by 2'^^^ and summed over j to obtain 



oo „ 



dv 



dv' 



d{v, v') 



l\n+2s 



{{v){v')) - ld(.,„o<i' (0,1). 
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The derivation of this inequahty is similar to that of Proposition [9] (which has the 
advantage of simpler notation), so we remind the reader of its proof now: 

Proof of Proposition\^ For any j > 1, one has the equality 

I fdvf dv' fdz ifiv) - f{v')fi,,{z-v)^,{z^v!_) {v) {v') {zY 

dv {\Q,f]{v)Y {vY ^ I dvf dz(f{v)fi^,{z-v)Q,(l){z){v){zY , 



simply by expanding the square (f{v) — fW))'^ ^-^^ exploiting the symmetry of 
the integral in v and v' (note also that the corresponding statement holds true for 
Q{) when ip is replaced by if). By our assumption on the projections Qj, namely 
|Qj(l)(z)| < 2~^^, we may control the second term on the right-hand side by 

2-2^' / dv {f{v)f{vY- 



This bound follows from the change of variables zi~>w + 2 '-^t^u, which is a change 
the variable from z to u and has Jacobian {v)~^ , so that 

dz\^jiz^v)\{zY ^{vy^ I du\'P{T^u + 2-^-^\u\^e,,+i)\{v + 2~^T,uY 

We have used that (v) « (z) on the support of \ipj{z_ — v)\ since j is non-negative. 
Moreover, the same change of variables can be used to show that 

dz \i^,{z-v)\\^,{z-v[)\ < 2"^((^'> W))"^- 

Here we use the inequality \ipj{z — i/)! < 2"^ and observe that (v) « (z) « {v') on 
the support of the original integral. Moreover, the triangle inequality guarantees 
that the integral is only nonzero when d{v,v') < 2^-'+^, so that we have 

2'^' f dz\ij,(z-v)\\i,,{z-vO\ {zY < {{v) {v'Y'i^ 22^^2"n,(„.,,)<2-.>i. 

These estimates may be summed over j > I (when v ^ v') because the sum 
terminates after some index jq with 2^-'" < d(v,v') < 2^^°+^, yielding ([5S|) . since 

oo jo 

^2^'*^2"-'lrf(i,_„/)<2-3+i = ^ 2^'*-'2"^ld(t,,t,')<2-3+i ^ d{v,v'Y"-^'^'''ld(vA,')<i- 
The last inequality follows from 2(^''+")-'" < (i(u, The remaining term 



j = is already bounded above by |/|?2. Q 



Thus, following the proof of Proposition [9l subject only to the establishment of 
the decay condition ([57]) similar to the hypothesis of Proposition [31 it follows that 

oo ^ 

(58) J2 E2'^^^'"'^'7 dv\dpy'^Q,fivr{vY^''w''-'\^\{v)<\f\%s., 

also holds for any multi-index of derivatives a on and any fixed K and £. The 

catch, of course, is that the functions tp'' become increasingly difBcult to control 
when either \a\ or K becomes large. This means that the uniform estimate is 
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increasingly difBcult to obtain. In the next section, we will establish the desired 
inequality contingent on the following cancellation condition; that 'ip^ satisfies 

dup(u)(V"V'')(z„u) = 0, 

for all polynomials p of degree 1 and all mult i- indices \a\ < 1. Since ■0*^ is itself 
related to the original ^ by taking a sequence of V-derivatives and multiplying by 
a polynomial, it is sufficient to choose the original so that ip satisfies 

(59) / d«p(«)(V"^)(r,«) = 0, 

for all polynomials of degree at most M and all \a\ < M, where M is any fixed but 
arbitrary natural number. 

3.3. Selection of (p and inequalities for smooth functions. Let D be the 



dilation on functions in M"+^ given by Dtp{w) tp . We choose a radial function 
tfo G C°°(M"+^) which is supported on the ball \w\ < R, with R > 0, and satisfies 

(60) / du ipoiZvU) = 1, Vw G M". 

Since ipo is radial this equality will be true for all v if it is true for any single v. If 
p is any homogeneous polynomial on R", a simple scaling argument shows that 

dup(«)(V"oi^^0)(T,7/)-2"+d'=SP-|a| / dup{u){\/"^o){T,u). 



Next, fix some large integer M and consider the C°° function ip which is supported 
on the ball of radius 2^^^ R and is given by 




By induction on M and the scaling argument just mentioned, it follows that (p 
satisfies exactly the same normalization condition as (/Sq, namely Assuming 
that M is fixed, the radius R may be chosen so that 2^^^ R < 1. For this fixed M, we 
take the corresponding function ip to be the basic building block of our Littlewood- 
Paley projections, i.e., we use this ip in the definition ([55]) . and use ip = ip—2~^Dip 
in ([M)) . By our particular choice of ip, this ip satisfies the cancellation conditions 
(|59)) for all polynomials p of degree at most M and all multi-indices a of order at 
most M . Consequently, we have the following integral estimates: 

Lemma 8. Choose a large integer M . Suppose that (|59p holds for all polynomials 
of degree at most M and all multi-indices a of order at most M . Fix any multi- 
index a with n + 1 components. If f is a smooth function on M", and g is some 
non-negative function satisfying for all v € M" that 

sup \idpfKv')\<Cf giv), 

d(«,"')<l|^|<fc+l 

Suppose k is some integer satisfying — 1 < fc < ^{M — \a\). Then the inequality 
(61) |2-I"I-'V"Q,/(«)| < 2-('^-+i)-''Q g{v), 

holds uniformly in f , v, and j > 0. 
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Proof. We first consider the case a = 0. We study the integral Ij{v) given by 
== / dv' ipj{y_— v^uj{v') 

= {vy'^ / du ^p{T^u + 2-^-^\ufen+i)ujiv + 2-^T^u), 

where the equality follows after the change w' i-> w + 2^^t^u. If we expand the 
integrand in powers of 2^-' by means of Taylor's theorem, we have an asymptotic 
expansion of this integral, with the coefficient of 2~''^ equaling 




This can subsequently be expanded as a sum of terms, each of which is an integral 
of some derivative of iJj times a polynomial in u. The order of differentiation is at 
most k, and the degree of the polynomial is at most 2k. Consequently, if 2k < AI, 
the fc-th term in the asymptotic series will vanish identically. Using the integral 
form of the remainder in Taylor's theorem, it follows that Ij{v) is exactly equal to 

(w)''-^ X J de {i-e)''-^[,p{T^u + 2-^-^9\u\\n+iMv + 2-^eTyu)], 

for any k < From here, it is elementary to see that 

|/.(«)l<2-('+'^M«>"' sup \^d0u:){v% 

uniformly for all v and all j > 0. Setting ijj{v') ~ {v') f{v') establishes the result 
for a ~ 0. When a ^ 0, the effect of including an additional 2~l"l-' V" is to replace 
ijj with V^ijj in the definition of Ij. The proof follows exactly as before, where now 
one only has ([5^ for derivatives up through order M — |a|. □ 

This lemma may now be directly applied to obtain the uniform bounds (|57p on 
\Qj{uji){v)\ as required for Proposition [9] and (f58| to be true. 

4. Upper bounds for the trilinear form 

In this section, we finish off the proofs of Lemma 121 Lemma 21 Lemma [5l and 
Lemma [SI from Section 11.51 The main tools we use in each of these proofs are the 
estimates from Section [2] with the decomposition of the singularity therein, and also 
a decomposition of the functions themselves via the non-isotropic Littlewood-Paley 
theory developed in Section [3] We begin with a variant of Lemma [5l 

4.1. The main non-linear estimates. We will now prove non- linear estimates 
in the velocity norms in Lemma |9l They are somewhat technical to state because 
we would like to distribute negative index decaying velocity weights among the 
post-collisional velocities ^ . This is explained above Proposition [31 We have 

Lemma 9. (Non-linear estimate) Consider the non-linear term ([9]) and (|22p . 

For any multi-index 13, any ,i~ > with i = i'^ — £^ and £' < £^ we have 

(62) \{w''r0ig,h)j)\< k'^-'5lLH/^U-,,J/U|- 

+ \w'^-''g\^.^Jw'+''h\r.^\fUs.^. 
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Here £ + £' ~ £^ — {£ — £'). We alternatively use the Sobolev embedding on g: 
(63) I {w^'Tp{g, h), f) I < \g\ ^. |/iU=.. 

+ \w'^-''g\^.^Jw'+''h\r^.\f\^^^^.. 

Of course, these estimates immediately imply Lemma[51 as we explain just now. 
Consider derivatives of the no n- linear term as in (j22p : a typical term is 

{w^'~'\^\Tp,{d';Z;^g,d;^^h),d^f) . 

Note that since \a\ + \/3\ < K we also have |ai| + |/3i| + |q! - ai| + |,3 - ^i| < K. 
To estimate this term, we will apply either (j62p or (|63p : in the former case, the 
right-hand side applies an additional K* = + Ij velocity derivatives to d^^h 
(a consequence of Sobolev embedding), and in the latter, the same K* velocity 
derivatives are applied instead to d^Zp^d- We will make the choice of ([5^ versus 
((63)) so that the total number of derivatives on either h or g does not exceed K . 

To that end, consider the situation when |q;i | + |/3i | < if — 2K^. Applying (p^ 
with 1+ = 1>Q, £' = |/3| and f = |/? - so - f = |/3i|, gives the estimate 

(w''-^\f'\Tp,{d'^pZtlg,d1lh),dp) 

Note that the £ right here is not the same £ as the one in Lemma O instead we 
use Lemma [HI with £ replaced by £ — Now after also integrating over T" we 
can use the Sobolev embedding 7/^"(T") C L°°(T") apphed to the norms of the h 
terms and a simple Cauchy-Schwartz to establish LcmmaO Here we have used this 
Sobolev embedding twice on the h terms, so the maximum order of differentiation 
of h on the upper bound wiU be |ai| + |/?i| + 2A'* < K. 

By symmetry, if ja — ai | + |/3 — /3i | < K ~ we may instead apply (|63p and 
estimate the g terms on the right-hand side via Sobolev embedding applied to the 
a;-variables as in the previous case. In the remaining situation, we have 

\ai\ + \l3i\>K-2KZ + l, 
\a-ai\ + \l3^/3i\>K-2KZ + l. 

However, since |ai| -I- |/3i| 4- jo; — aij 4- |/3 — /3i| < K, it must also be the case that 
\ai I + 1^1 1 < 2K* - 1 and |a - ai I 1,3 - ,3i I < 2K* - 1. In this situation we wiU use 
([5^ (the Sobolev embedding for the velocity variables is applied to h) and estimate 
the g terms on the right-hand side via Sobolev embedding in the x- variables. For 
both g and h, the total number of derivatives appearing on either term will then 
be at most (2K* — 1) -I- K*. Therefore, Lemma[5]will hold whenever K > 3K* — 1. 
Now we set about to prove the two non- linear estimates in Lemma [9l 

Proof of Lemma\^ We'll write 

oo oo 
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and likewise for h. Then we expand the non-hnear term as in ([22]) as fohows: 

CXD 

{rp{g,h)j) = ^ (r^(5,v),/,) 

hj'=0 

OO oo oo 

(64) = J2 (r^ff: h,), /,) + ^ 5] {(r^(5, h,+i), /,) + {T^{g, h,), . 

j=0 1=1 j=0 

Consider the sum over I of the terms {Tp{g,hj^i), fj) for fixed j. We expand T 
by introducing the cutoff around the singularity of b in terms of tJJ'^ and t'I'^: 



oo OO 



^{^p{9,h.j+i)Jj) = ^ ^^Tl'\g,h.j+ufj)-T'l-\g,hj+i,fj)^ 

1=1 k=-oo 1=1 



(65) = J2 {Tt^\9,h~PjhJ,)-T'i^'(g,h~P,hJ,)} 

k— — oo 

oo oo 

(66) {^+'(9, h,+i,.U) - Tl'^ig, h,+ij,)} . 

1=1 k=l 

Here we have used the basic telescoping property h — Pjh ~ T^^i^j+i- It is 
worth remarking that throughout the manipulation, the order of summation may 
be rearranged with impunity since the estimates we will employ below imply that 
the sum is absolutely convergent when g, h, f are all Schwartz functions. 

We will prove first , and later explain how to similarly obtain ([5^ . Regarding 
the terms the inequalities and dictate that 

\Tl''ig, h - P,h, /,) - Tl'\g, h P,h, f,)\ 

fc— — OO 

Notice that the upper bound in ([5^ is bigger than the upper bound from . We 
have used \w'^ Pjh\i^2 < \wPh\]^2 Vp G M, a consequence of ([55)) . Thus 



Y \Tl-'ig,h-P,hJ,)^Tl-'ig,h-P,h,m 

j=0 fc=— OO 



< 



i=0 



< 



This is just Cauchy-Schwartz. The favorable comparison of the square-function 
norm of / to the norm I/Iat^ ^ is provided by ([5^ 
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As for the terms (|66|) . when k < j a similar approach holds; namely, ([24]) and 
35| guarantee that 



fc=i 



|r^'^(5, hj+ijj) - T^'\g, hj+ijj) 



(we have used the trivial facts that 2^'''^ = 2^''i2^'''-''-^^ and Y^l^^ 2^"''^-^^ < 1). In 
particular, this inequality may be summed over j; another application of Cauchy- 
Schwartz gives that 



j=0 



1+ -I' 



3=0 



2 

J\Ll 



3=0 



This estimate may clearly also be summed over Z > 0. 

A completely analogous argument may be used to expand T p for the terms in 
(jM]! of the form {Tp{g, hj), fj+i) in terms of T^''^ - T^^; 

OO OO OO 



l — l k — — QO l—l 

(67) = {Tl-'-Ty){g,h,J~P,f) 

OO OO 



1=1 k=l 



In this case the estimates and ([5^ are used to handle the terms (|67p just as 
the corresponding terms (|65p were handled. Regarding the sum now and 
([55]) are used to estimate the sum analogously to the estimates of for k < j. 
The only difference is that the roles of h and / are now reversed. 

Recalling the original expansion of IJ^ fj{g^h), f) it is clear that the only terms 
that remain to be considered are the following: 



OO CJO OO 



(69) f,{g,hJ)''^Y.Yl E {Tl-\9.h,+uf,)-Tl-\g.h,+ufj)} 

1^0 j^O k^j + 1 

OO OO OO 

(70) +EE E {Tl'\g,h,J,+i)-Ty{g,h,J,+i)]. 



1=1 j=0 k=] + l 
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In other words, wc have ah'eady estabhshed the inequahty 



{Tpig,h)J)-T^ig,hJ) 



< 



e+e' , 



I/I 



Now the terms on the right-hand side of and the terms ([Tn]) are both treated 
by the cancehation inequahties. The terms (|69)) . for example, arc handled by (|4T|) . 
For any fixed l,j,k, wc have 



There is decay of the norm as A: — oo since 2s — 2 < 0. This yields 



k=j+i 



^(5, hj+ijj) ~T'l'\g, h-i+ijj) 



< 2(2«-2)i| 



'H 



Iw'h 



Again Cauchy-Schwartz is applied to the sum over j. In this case 2'^^'' ^^-^ is written 
as 2(''-i)-''2''(-'+')2-'''2--''; the first factor goes with fj, the second with h^+i, the 
third and fourth remain for the sums over I and j respectively. Once again ([55)1 is 
employed. The factors of 2"*' and 2^-' allow us to finish the sums over I and j. 

The desired bound for the non-linear term is completed by performing summa- 
tion of the terms ([70|) . The pattern of inequalities is exactly the same as the one 
just described, this time using (|46| . In particular, one has that 



Tl\g,h,J,+i) -Ty{g,h,J,+i)\ 

<2(^-^)^|5l^_.r.kV.+/lL^, 



iVK/i, 



Again, since 2s — 2 < 0, this leads to the following inequality for the sum over k: 



k=j+i 



\T+\g: hj+ij-i) ~T'l-\g, hj+ijj) 



The same Cauchy-Schwartz estimate is used for the sum over j; there is exponential 
decay allowing the sum over I to be controlled. In the estimates above m > can 
be taken arbitrarily large, yet — m — iis sufficient. The end result is precisely (|63l) . 

We can obtain using the same techniques. We expand all of the sums from 
((64)) in the same way. Then we estimate (|65p using the inequalities (|25l) and ([33|l . 
We control ([Ml), when fc < j, with ^ and Then jSO]) and ^ are used to 

handle the terms ([67]) . Now ([30]) and (|36|) are used to estimate the sum ([68]) for 
k < j. For the cancellations, (|69| arc handled by (|40|) . Lastly, (|70|) is controlled 
using (|45p . Otherwise, the pattern of summing these inequalities is exactly the 
same as the one just described above. □ 



This concludes our main non-linear estimates. 
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4.2. The Compact Estimates. Here we collect some estimates for the linearized 
collision operator. The first one is the key to the estimate for operator K in ([T5|) . 

Lemma 10. (Compact Estimate) For any ^ e R, we have the uniform estimate 

\{w^'Kg,h)\ < \w'g\,^2^^^_yh\L2^^^_^, 6 ^ mm{2s, {n - 1)}. 

Since (5 > above Lemma [TUl easily implies ([T5|) when g ^ h. To see this, first 
apply Cauchy's inequality with ^ to the upper bound in Lemma [TOl 

\{y^''Kg,g)\<^\w'g\l2 +C\-,\w'g\l. 

For the term C„\w^g\'j-2 above, we split into \v\ > R and \v\ < R. Choosing 

^■, + 2s^S 

R > sufficiently large so that C,jR^^ < ^ proves ([T5|) subject only to Lemma [TOl 
Now Lemma [10] and other estimates will follow from (f72|) below. 

Proposition 10. For any function (j) satisfying (j23p . we have the estimate 

(71) \{w''rpi^,h)j)\<\h\^s..\f\^s... 

For the next two estimates, we suppose that (j) further satisfies (|47p . Then 

(72) |(^2%(.g,</>),.f)| < h'5lL^,,,_^„_Ju;VlL^^,^_,„_,,. 
For any m > we also have 

(73) |(z«2^r^(.9,/),0)|<|.9|i2_|/|i2_. 

Each of these estimates hold for any /3, and any ^ € R. 

Notice that these imply several other previously-stated estimates. In particular, 
Lemma [TUl is a special case of ([7^ and Pao's estimate of i^Kiv) in ([T^ . Thus 
Proposition [TOl implies Lemma [3l since also ([T4| follows directly from (j7ip . Other 
uses of Proposition [To] will be seen below. 

Proof of Proposition] in To prove ([TT]) . we expand (w^^r/3(0, h), /) as in ([64]). and 
the proof follows the same lines as the proof of ([55)1 . Following that proof, we 
estimate ([65]) and ([66)) . when k < j, using the inequalities ([24)) and ([52]) . Then ([29]) 
and ([52]) arc used to handle the terms (|67p and (|68p for k < j. For the cancellations, 
([69| are handled by ([4T|) and ([70| is controlled using ([46]) . The pattern of summing 
these inequalities is exactly the same as the proof of ([63| . 

To prove the estimate in (|72p we will use the inequality (j6ip . In particular 



(^2%(g,0),/)| 



E 

j^O A;— — oo 



oo 

^ {T^^\g,0„/)-T^^^(3,0„/)} 



^ l-'^^l^?+.-<„- J-Vl.^,_,„_, E E ^fin{2(--)^2-}2-^■. 

J— k— — oo 

We have used ([26]) and ([50]) with Proposition [7] Specifically, in each of those 
estimates 4)j satisfies ([47]) with < 2^^-', as in ([6T]) . 

The estimate for (uP'^T p{g, /), 0) in (|73p is proved in exactly the same way using 
instead ([26]) . ([ST]) . and Proposition [7] In particular. Proposition [TOl follows. □ 

This concludes our compact estimates. 
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4.3. The coercive estimates. In this sub-section we will prove the coercive in- 
terpolation inequalities in (fTB]) and (fT7|) from Lemma El and Lemma ID Actu- 
ally, P7)) is a trivial consequence of Lemma U] and because {w'^^Lf, /) = 
{w^^Nf, f) + {w^^Kf, /). Thus we will restrict attention to p6)) and Lemma S] 

Proof of Lemma Firstly, Lemma |4] for £ = and C = is proven in [331 Lemma 
4] ; note that |33l Lemma 4] is only stated for n = 3 dimensions and a more restrictive 
class of collision kernels than ^ and ([3]), however the proof directly extends to the 
case of n > 2 dimensions and all of the collision kernels 7 -|- 2s > — (n — 2). We 
focus here on estimating the norm piece when £ ^ 0. We expand 

{w''Ng,g) = \g\% + f dv w^\v)u{v) \g{v)\^ + J, 

where \g\Bi is defined in (|12p . Furthermore, we have the following equivalence 

\9\\ + Wg\l.^^^^\g\l^,-.. 

The lower bound > of this equivalence follows directly from the proof in [331 Section 
6]. The upper bound, <, follows from the estimates for (w'^'T{M^ g), g) in ([71]) when 
combined with the interpolation inequality that we prove just now. 
The error term, J, above takes the form 

J=\ I dv f dv, [ daB{g' -g)g{w''\v')^w''\v))M'JU. 

These expressions are derived exactly as in the computations preceding (1111) . 

We will show that this error term J is lower order via an expansion of the kernel. 
In particular we claim that there exists an e > such that 

\J\<\9Ww'g\L2^ . 

This estimate implies that J satisfies the same upper bound as ([T5|). which easily 
implies Lemma [H by following the procedure which is explained below Lemma 1101 
To prove this claim, notice that Cauchy-Schwartz gives us 

/ /■ r f (w^'^(v') -w^'^(v))^ V^^ 
\J\ < \9\b, / dv dvJ daB\g\'^ ^ , ^ M'^M, . 

We will show in particular that there exists an e > such that 

/7A\ f ^ f ^ U( {w'^'^{v')-w'^'^{v)f 2£/ X/ \7+2s-. 

(74) / dv, I da B{v-v,,a) jj-- M,M,<w [v){v)' 

and this will establish the claim. 

To obtain ([74p . first a simple Taylor expansion yields 

w^\v') - w^\v) = (v' - v) ■ (Vw2^)(7(t)), 3t e [0, 1], 

where 7(t) = v + t{v' — v). Since \v' ~ v\ = \v', — v,\ we have 

(7(r)) < (v) {v' ~v)< (v) « - V,) , 

and similarly, {j{t))~ < (v)" {v', — v,) . Thus generally 
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Now we split \v' - = \v' - - v^]'^-^"-^ for any (5 G (0, 2 - 2s). We can 

expand \v' — v\'^''^^ = b ^ w*^'^'' (sin|)^''^'' , and then we clearly have 

/ da B{v~v,,a) \v'~if'+^ < \v - v,\-y+'''+^ . 

This follows directly from ^ and (g]). Furthermore, |< -i;*p"2s-5(-^,j/^,^^)i/4 < ^ 
Putting all of this together, we see that (f74| holds with e ~ 2 — 5 > 0. □ 

With the help of our coercive estimate ([T7| with no derivatives, in the following 
we will prove the main coercive estimate with high derivatives. 

Proof of ([T6| . We use the formula for Lg from ([T0|) . As in p2p . we expand 

After multiplying by w'^^~^^^^d"g, and integrating over R" we can estimate the term 
^y^2£-2|/3|^ (^d^g^ ^d^g) as in (fT7|) . For the error term that arises, i.e., \dpg\L'^(Bii) 
for some i? > we use the compact interpolation for any small 5 > Q: 

\dpg\L^{BR) < v\dpg\HHB„} + C\9°'9\lhBc)- 
Here C > is some large constant, and jy > is any small number. Further 
\9^9\h^{Bp{) ^ \9pg\N^-''^^^ when S < s; this holds because the non-isotropy of the 
norm N^'^ only comes into play near infinity. More precisely, if v and v' are 
confined to the Euclidean ball of radius R at the origin, then d{v,v') ^ \v ~ v'\ 
(with constants depending on i?), and so on this region the expression for ([7]) is 
comparable to the Gagliardo-typc semi- norm for the space H^{Bft). This gives the 

estimate for L {d'^g^ ■ 

We estimate the rest of the terms using (|7ip and ([7^ . We use the extra velocity 
decay which is left over from these estimates to split into a large unbounded region 
times a small constant and a bounded region with a large constant. On the bounded 
region, we use the compact interpolation as just used in the last case to put all of the 
velocity derivatives into slightly larger Sobolev norm multiplied by an arbitrarily 
small constant. This is all that is needed to finish the estimate. □ 

5. De-coupled space-time estimates, global existence and rapid decay 

We will now proceed to show that the estimates proven in all of the previous 
sections can be used to establish global existence as stated in Theorem [TJ The 
methodology that we employ, which goes back to Guo [37], essentially de-couples 
the required space-time estimates that are needed from the new fractional and non- 
isotropic derivative estimates which are shown in the previous sections. This works 
precisely because of the specific structure of the interactions between the velocity 
variables and the space-time variables. Our key new contribution in this direction 
is to identify the unique weighted geometric fractional Sobolev norm ([7]), and prove 
sharp non-isotropic geometric fractional derivative estimates such as those stated 
in Lemma m Lemma [U Lemma El and Lemma El These can then be included into 
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the general de-coupled energy method in order to prove global existence and decay. 
We point the reader's attention also to the general abstract framework of [61]. 

Wc furthermore remark that the whole space case of our main Theorem [1] when 
X S M" replaces x € T", will follow from our new velocity estimates in the previous 
sections when combined with other known whole space cut-off methods. This would 
only require re-writing this last section, using instead the cut-off energy methods 
in the whole space such as for instance j38 l l41 ] [52] . 

We will initially discuss the coercivity of the linearized collision operator, L. 
With the nuU space and the projection we decompose f{t, x, v) as 

(75) / = P/ + {I - P}/. 

We prove a sharp constructive lower bound for the linearized collision operator. 

Theorem 11. There is a constructive constant > such that 

{Lh,h) > 6o\{I~P}h\%.,,. 

This coercive lower bound is proven with our new constructive compact estimates 
from Lemma [3] and Lemma [4] when used in conjunction with the non-sharp but 
constructive bound from Mouhot [48j for the non-derivative part of the norm. 

Proof of Theorem llll Suppose /i = {I — P}h. From p7p . for some small ?; > 

{Lh, h) > ^^\h\l.^., ~ C\h\l,(^B^,y 3C > 0. 

Here Be is the ball of radius C centered at the origin. The positive constant C is 
explicitly computable. From [48], it is known that under our assumptions 

{Lh,h)>5i\h\l,. 

Here (5i > is an explicitly computable constant, and 7 is from (|5]). 
Lastly, for any (5 G (0, 1) we employ the splitting 

{Lh, h) = S{Lh, h) + {l- S){Lh, h) > 57j\h\%,., - + (1 - 5)5i\h\l2. 

Since C is finite notice that \h\'j^2(^Bc) ^ l^lt^ 7 G M. Thus the lemma 

follows by choosing 6 > sufficiently small so that the last two terms are > 0. □ 

5.1. Local Existence. In this section, we explain the approach from our first 
paper |33| to prove the a priori energy estimates which are the key step in the local 
existence theorem. Our local existence proof for ([8]) with small data is based on a 
uniform energy estimate for an iterated sequence of approximate solutions. 

Our iteration starts at x, v) =' fo{x, v). Wc solve for f^'^^{t, x, v) such that 

(76) (9t + i;-V,+iV)/'=+i+A7''^=r(/\/''^+i), f''+^{Q,x,v)=fo{x,v). 

It is standard to show the linear equation (|76|) admits smooth solutions with the 
same regularity as a given smooth small initial data, and also has a gain 
of L'^{{0,T); Ng''^). This does not create difficulties and can be proven with our 
estimates. Wc explain herein how to establish the a priori estimates necessary to 
ffnd a local classical solution in the limit as fc — > 00. The key point will be to obtain 
a uniform estimate for the iteration on a small time interval. The crucial energy 
estimate is as follows: 
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Lemma 12. The sequence {f''{t,x,v)}k>o is well-defined. There exists a short 
time T* > 0, such that for \\fo\\H sufficiently small, there is a uniform constant 
Co > for which 

sup sup f||/^t)||l,+ fdr Wfirnl) <2Co\\forH- 

k>0 0<t<T* \ Jo J 

The proof of this estimate is exactly the same as the proof in [33] using the 
estimates in Section [TTSl see [33] for the details. With our uniform control over the 
iteration ([76]) granted by Lemma|T2l we have the following local existence theorem. 

Theorem 13. (Local Existence) For any sufficiently small A/q > 0, there exists a 
time T* ~ T*(Mq) > and an Mi > 0, such that if 

II/oIIh<a/i, 

then there is a unique solution f(t,x,v) to ^ on [0,T*) x T" x M" such that 

sup g{f{t))= sup (\\fml+ f dT\\f{T)\\%) <Mo. 
0<t<T' 0<t<T' \ Jo / 

Furthermore is continuous over [0,T*). Lastly, we have positivity in the 

sense that if Fo{x, v) = /i + i-J-^^^ fo > 0, then F(t, x, w) = /i + i-i^^^ f{t, x, v) > 0. 

For the proof we refer to the argument in [33], the case here is the same. 



5.2. Coercivity estimates for solutions to the non-linear equation. The 

next step is to prove a general statement of the linearized H-Theorem, which man- 
ifests itself as a coercive inequality. These types of coercive estimates for the 
linearized collision operator were originally proved by Guo |36[l37j in the hard- 
sphere and cut-off regime. In [33], we proved this estimate for the Boltzmann 
equation without any angular cut-off, as in ([3]), but with p > 3 or 7 -f 2s > 0. The 
next theorem extends this estimate to the full range of p > 2 and more generally 
7 + 2s > —(71 — 2) in 71 dimensions. 

Theorem 14. Given the initial data /q £ Hf (T" x M") with derivatives of order 
K > K(n) and weights of order £ > 0. Suppose that fo satisfies ([6]) initially and 
the assumptions of Theorem \13[ For the corresponding local solution, f{t,x,v), to 
([8]) which satisfies ([6|, there exists a small constant Mo > such that if 

(77) ||/W|Ph<Mo, 

then, further, there are constructive constants S > and C2 > such that 

\a\<K \a\<K 

where T{t) is the "interaction functional" defined precisely in (|80p below. 

This Theorem [14] is proven in exactly the same way as in |33| , using the usual 
decomposition of P/'s coefficients a-f{t, x), b{{t, x) and c-f{t, x), for a solution to (jS]) 
using (Unj). We explain briefly how this works and refer to [53] for the full details. 
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We can decompose a solution / as ([75]) to obtain the macroscopic equations 

dtc + ^ib^ = -dtTi + k + Ti 
dihj + djh = -dtrij + kj + T^- {i ^ j) 
dth + diQ = -dtm + hi + ^bi 
dta = -dtra + la+T^a- 
Above di = dxi- Consider the basis, {e^;}^,"'!^^'^"''""^''^^, which consists of 

(78) {vf^)^^,^^^^, {v^vj^/r^)i<,<j<r., («»Vm)i<,<„, Vm- 
Further, for notational convenience we define the index set to be 

M = |c, i, , bi, a| i, j = 1,. . . 

This set A4 is just the collection of all indices in the macroscopic equations. Then 
for 1^1 E Ai we have that the l^{t,x) are the coefhcients of — P}/), given by 

(79) P}/) -V ■ V,({I - P}/) - L{{1 P}/), 

with respect to the elements of (j78p : similarly, the r^ (i,a:) and r^ (t, x) are the 
coefficients of r(/, /) and {I— P}/ respectively. Each r^, for example, can be 
expressed as 

3n+H-n(ri-l)/2 

E cn{i-p}/,e,). 

fc=i 

All of the constants can be computed explicitly. Furthermore, each of the terms 
and can also be computed similarly. 
Now, wc define the total interaction functional as 

(80) x{t)'^ J2 {i:{t)+ib{t)+icit)}, 

\a\<K-l 

which is composed of the following elements: 

X;^{t) = f dx {W- <9"fe) d°'a{t, .t) + V / dx d^d'^n, d°'a{t, x), 

I^it) = - / dx d"r,{t, x) ■ V,a"c(i, x). 

This interaction functional will be comparable to the energy defined below. 

With these macroscopic equations, the key to proving Theorem [14] will be the 
following lemma. In fact, once we have this next Lemma [T5| the proof of Theorem 
[THwill follow exactly as in [33] and we refer to that paper for the rest of the details. 

Lemma 15. For all of the microscopic terms, Z^, from the macroscopic equations 

Ell^A^llfff-^ E !I{i-p}9"/||l^^,^(T"xr")- 

MeA4 \oi\<K 
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Furthermore, let {71^ he valid for some Mq > 0. Then for K > 2K* we have 

t^eM \a\<K 

We are using the abbreviated notation Ll = L2(M") and iJf = iJ^'(T"). 

Proof. Recall {ek}, the basis in (|78p . To estimate the /^j terms, it suffices to estimate 
the iJf -1 norm of (/({I - P}/), e^). We use ^ to expand out 

- P}/), efc) = -(^; • V,({I - P}a"/), e,.) - - P}5"/), e^). 

Now for any |a| < ^ ^ 1 by Cauchy-Schwartz we have 

\\{vV,i{l-P}dy),ek)\\h < I dxdv \ek{v)\ \v\^ |{I - P}V,a"/|2 

<||{I-P}V,5"/||i.(T"xR")- 

Here we have used the exponential decay of efc(w). 

It remains to estimate the linear operator L. Notice that the basis vectors Ck 
satisfy (|47)) . With the expression from (|10p and the estimate (|73l) . we have 

\{Li{I-P}dV),ek)\\l.< |{I-P}a"/|^. |M| 



L2 



This completes the proof of our estimates for the l^ terms. 

As in the above, to prove the estimate for F^, it will suffice to estimate the 
norm of (r(/, /), efe). We apply ([75)) from Proposition [TUl to see that for any m > 

|q|</<' C(i<a 

We take the supremum over the term with fewer derivatives to obtain for K > 2K^ 

\a\<K 
\a\<K 

The last inequalities follow from the Sobolev embedding L°°(T") D iy''^"(T"). □ 

We are now ready to prove global in time solutions to (|5]) exist. 

5.3. Global Existence and Rapid Decay. With the coercivity estimate for non- 
linear local solutions from Theorem [TH we show that these solutions must be global 
with the standard continuity argument which is used in [33]. Then we will prove 
rapid decay for the soft potentials, using the method from the cut-off case [53]. 
A crucial step in this analysis is to prove the following energy inequalities: 

(81) + T^i,mit) < Ce^rnyMU)Ve{t). 

These hold for any £ > and m = 0, 1, . . . , A'. We define the "dissipation rate" as 



\/3\<m \a\<K-\l3\ 
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We also write V^.k — Vi ^ \\f{t)\W- Furthermore the "mstant energy functional" 
£i,m(t) for a solution is a high-order norm which satisfies 

|/3|<m \a\<K-\p\ 

Wc prove this energy inequality for a local solution via a simultaneous induction 
on both the order of the weights £ and on the number of velocity derivatives m. 

The first inductive step is to prove (|8ip for arbitrary spatial derivatives with 
£ = and |/3| = 0. We first fix Mq < 1 such that both Theorems [HI and [Til are 
valid. We now take the spatial derivatives of d" of ^ to obtain 

(82) E (^9"/,a"/)= J2 id"rifj),dv). 

\a\<K \a\<K 

By Lemma Owe have 

^ (5"r(/, /), a"/) < y/s^ Vo{t). 

\a\<K 

Now with Theorem [TT] and then Theorem [14] we have 

^ (La"/,a"/)>5o ^ i|{i-p}a"/i|| 

la|<i<: |q|<A' 



2 

ATS, 7 



2 



|a|<X |a|<A' 

With ^ min {f , > and C" ""^^ 5qC2 > 0, we conclude that 

Now, by (j80|) . for any C" > we can choose a large constant Ci > such that 

WfmiiH^ < (Ci + 1) wfmiiH^' cm < ii/wiii.^.. 

Notice Ci only depends upon the structure of the interaction functional and C", 
but not on f(t,x,v). We then define the equivalent instant energy functional by 

£oAt) = iCi + i)\\.fmliH'<-c'i{t). 

We multiply ([5^ by Ci and add it to this differential inequality to conclude 

+ ~ST^oAt) < a v^Po(t), a > 0. 

at 

In the last step we have used the positivity of L > 0. We have thus established 
(t8T|) when £^\I3\^ 0. 

We turn to the case when £ > 0, but still \(3\ = 0. We only have pure spatial 
derivatives. With pT)) in Lemma [6l we deduce that for a C > and £ > 

(83) {w''LdV,dV) > - 

Take the derivative of ([S]), then take the inner product of both sides with w^^d"f 
and integrate to obtain the following: 

'ld_ 
2 dt 



E ( ^T^d-fmh + {w''Ldv,d-f) ) < 



y\<K 
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We have used Lemma [5] to estimate the non-hnear term. We apply the coercive 
lower bound (|83|) . Then we add (|8T|) for the case £ = |/3| = multiplied by a 
suitably large constant C2 to the result. This yields 

where £e.o{t) = \ '}2\a\<K Il^^^"./(^)lli2 +^^2^^0,0(0- Since this is indeed an instant 
energy functional, we have (|8ip when ^ > and — 0. 

The final step is the case when ^ > 0, but also |/3| = m + 1 > 0. We suppose 
that ([HI]) holds for any ^ > and any |^| < m. We take 5^ of © to obtain 

(84) {dt^v V,) dp + dpLd^f - - ^ < dp,v ■ y,di_pj + 9^r(/, /). 

l/9i| = l 

We use Cauchy's inequality for > 0, since = 1 we have 

Estimates using this particular trick were already seen in |36| . 

Now wc multiply (|84l) with u;^^~^l^l9"/ and integrate. We estimate the non- 
linear term of the result with Lemma [5j With (|16p we estimate from below the 

linear term inequalities above, we conclude 

([5T|) for 1/3 1 = m + 1, but only after adding to the inequality a suitably large 
constant times ([5T|) for \/3\ = m similar to the previous cases. This establishes 
(|8ip in general by induction. From here we can conclude global existence using the 
standard continuity argument that was described in Section 7.3 of jSSj . It remains 
to establish the time decay rates for the soft potentials, using the argument from 
[55] . Exponential time decay for the soft potentials, as in [SI], may also be feasible. 
If |l/n|l|f is sufficiently small, from (|8ip for m = K and ^ > 0, we have 

^SfXt) + 5Vi{t) < 0, 3,5 > 0. 
dt 

For the hard potentials, 7 + 2s > 0, exponential time decay follows directly from 
T^e.{i) ^ But for the soft potentials 7 + 2s < 0, the problem is that for fixed 

£, the non-derivative part of the dissipation rate 'Dg{t) is clearly weaker than the 
instant energy £e{t). In particular, we only have T>e{t) > £g-i{t). 

We interpolate with stronger norms to overcome this difficulty. Fix £ > and 
m > 0. Interpolation between the weight functions and yields 

The last inequality follows from Si^mit) ^ ^i+mi^)- Then for some Ci^m > 0, 

It follows that -m d{£i,{t))~'^/'^ / dt < -Ce,,n {£e+„i{0)y^^"' . Integrate over [0,t]: 
m(£-,(0))-^/'" - m(£-,(t))-^/" < - (0))"'/" C,,™ t. 
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Hence 

m))-'^"' > (f.+™(o))-^/™ + {f£(o)}-i/™. 

m 

Since we can assume ^^(O) < £g+m{0), the rapid decay thus fohows. Q.E.D. 
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